JID-A182  594  INFRARED  SENSOR  DESIGN  FOR  GEOSVNCHRONOUS  ANTISATELLITE 
VEHICLES(U)  AIR  FORCE  INST  OF  TECH  URIGHT-PATTERSON  AFB 
OH  SCHOOL  OF  ENGINEERING  J  G  HORNE  DEC  84 
UNCLASSIFIED  AF I T//GS0/PH/84D-4  F/G  22/1 


1/2 

ML 


lwrKAKtaj  SENSOR  DESIGN  FOR 


GEOSYNCHRONOUS  ANTISATELLITE  VEHICLES 


THESIS 

Jeffery  G.  Horne 
Captain,  USAF 

AFIT/GSO/PH/84D-4 


;-x  JUL2  7  w 

DEPARTMENT  OF  THE  AIR  FORCE 

AIR  UNIVERSITY  ^ 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


_  Wright-Patterson  Air  Force  Base,  Ohio 

i£-3  has  baea  uppromd 

ics  end  sal#;  fa  7  O  V  O 

'iic-triScScn  la  anUm»edL,.  U  ^  »  **  ^ 


AFIT/GSO/PH/84D-4 


INFRARED  SENSOR  DESIGN  FOR 
GEOSYNCHRONOUS  ANTISATELLITE  VEHICLES 


THESIS 

Jeffery  G.  Horne 
Captain,  USAF 

AFIT/GSO/PH/84D-4 


Approved  for  public  release;  distribution  unlimited 


AFIT/GSO/PH/84D-4 


INFRARED  SENSOR  DESIGN  FOR  GEOSYNCHRONOUS 
ANTISATELLITE  VEHICLES 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 

/ 

Requirements  for  the  Degree  of  v  \ 

>  % 

Master  of  Science  in  Space  Operations 


Accession  For 

NT  IS  GF.UI 
DTIC  TAB 
Unar.ic  ..  eed 
Ju  U.il'  Leation 


Jeffery  G.  Horne,  B .  A. ,  M.S 
Captain,  USAF 


By. 


Distri 

Avail 

Dist  I 


but  ion/  _ 

utility  Codes 
Av ail  and/or 
Special 


December  1984 


Approved  for  public  release;  distribution  unlimited 


In  performing  the  investigation,  analysis  and  writing 
this  thesis,  I  have  had  a  great  deal  of  help  from  others. 

I  am  much  indebted  to  my  faculty  advisor,  Major  James  J. 
Lange,  for  his  ideas,  suggestions,  directions,  and 
corrections.  I  would  also  like  to  thank  my  wife  Valerie, 
for  her  much  needed  support  and  assistance. 


Page 

ii 


Acknowledgments 

List  of  Figures .  v 

List  of  Tables .  vi 

Abstract .  vii 

I.  Introduction  .  1 

II.  Background  .  4 

Infrared  Sky  Surveys  .  4 

Infrared  Space  Sensors  .  7 

III.  The  Infrared  Sky .  12 

Walker  Sky  Model  . .  13 

Model  Data .  19 

Model  Results .  23 

IV.  Sensor  Design .  25 

Design  factors  .  28 

Target  Current  .  35 

Background  Current  .  35 

Signal-to-Noise  Ratio  .  36 

Maximum  Range  Estimates  .  43 

V.  Model  Accuracy  . .  44 

Limiting  Irradiance  Effects  .  44 

Comparison  to  Cataloged  Observations  .  .  48 

VI.  Additions  to  the  Model  .  . .  57 

Asteroids .  57 

Zodiacal  Light  Particles  .  61 

Revised  Maximum  Ranges  .  67 

VII.  Conclusions  and  Recommendations  .  71 

Appendix  A:  Coordinate  Systems  and  Conversions  73 

Appendix  B:  Cataloged  Observations  .  80 


iii 


Flux  Units 


92 


Appendix 


C: 


Appendix  D:  Walker  Sky  Model  Derivation  .  .  .  . 


94 


Bibliography  .  99 

Vita .  103 


List  at  Figures 


Figure  Page 

1.  Model  of  the  Galaxy .  14 

2.  Radiance  Contour  Map .  24 

3.  Sensor  Configuration  .  29 

4.  Sensor  Reticle .  29 

5.  Resolution  Elements  .  33 

6.  Number  of  Sources  versus  Irradiance  ....  46 

7.  Change  in  Radiance  versus  Irradiance  Range  46 

8.  Number  of  Sources  versus  Irradiance  ....  47 

9.  Change  in  Radiance  versus  Irradiance  Range  47 

10.  Selected  Sky  Areas  in  Galactic  Coordinates  50 

11.  Number  of  Sources  per  Square  Degree  ....  54 

12.  Ecliptic  in  Galactic  Coordinates  .  59 

13.  Radiance  versus  Elongation  Angle  .  64 

14.  Range  versus  Background  Radiance  .  68 

15.  Revised  Radiance  Contour  Map  .  70 

A-l.  Equatorial  Coordinates  .  78 

A-2.  Equatorial  Coordinates,  Horizon  View  ...  78 

A-3.  Ecliptic  Coordinates  .  79 

A-4.  Galactic  Coordinates  .  79 

D-l.  One  by  One  Degree  Volume .  98 


Table 

List  of  Tables 

Eag& 

I. 

Constants  for  Sky  Model  . 

17 

II. 

Number  of  Sources  Per  Square  Degree  .... 

19 

III. 

Total  Radiance  . 

23 

IV. 

Sensor  Parameters  . 

34 

» 

1 

V. 

Background  Current  . 

41 

VI. 

Maximum  Ranges  . 

43 

i 

VII. 

Centers  of  Selected  Sky  Areas  . 

49 

VIII. 

Radiance  Comparison  . 

52 

IX. 

Sources  per  Square  Degree  Comparison  .  .  . 

52 

1 

X. 

Number  of  Asteroids  per  Square  Degree  .  .  . 

60 

XI. 

Radiance  from  Asteroids  . 

61 

• 

XII. 

Zodiacal  Emission  1  . 

63 

XIII. 

Zodiacal  Emission  2  . 

63 

XIV. 

Revised  Maximum  Ranges  . 

67 

C-I. 

Magnitude  Difference  .  ..... 

93 

1 

vi 

l:'"  '  " 

AF I T/ GSO/ PH/ 8  4D- 4 


Abstract 

.  \ 

The  maximum  distance  from  which  a  possible 
antisatellite  sensor  can  reliably  detect  its  target  was 
investigated  based  upon  the  celestial  background  radiation 
at  a  wavelength  of  11  microns.  A  sensor  design  was 
postulated  and  assessed  with  various  background  conditions. 

The  analysis  was  first  conducted  using  the  Walker  Sky 
Model  to  estimate  the  amount  of  infrared  background 
radiance  from  stellar  sources.  Then  this  estimate  was 
compared  with  cataloged  observations  for  the  Air  Force 
Geophysics  Lab  catalog  and  preliminary  data  from  the 
Infrared  Astronomical  Telescope.  In  addition,  the  amount 
of  infrared  radiation  from  asteroids  and  zodiacal  light 
particles  was  also  considered. 

The  results  of  the  analysis  indicates  that  for  a 
sensor  using  a  single  detector,  a  spinning  reticle,  and 
optics  20  centimeters  in  diameter,  the  maximum  range  that  a 
one  square  meter  target  can  be  detected  reliably  is  between 
1300  and  2300  kilometers.  However,  more  investigations 
into  the  infrared  stellar  and  zodiacal  light  background  are 
needed  as  the  currently  available  data  is  incomplete. 
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INFRARED  SENSOR  DESIGN  FOR  GEOSYNCHRONOUS 


ANTISATELLITE  VEHICLES 

i.  intigdugtifln 

In  order  for  an  antisatellite  (ASAT)  vehicle  to  be 
effective  against  a  satellite  target,  it  must  first  detect 
the  target.  This  means  detecting  the  target  at  a  large 
enough  distance  away  to  allow  the  ASAT  to  maneuver  into  an 
orbit  that  will  result  in  interception.  The  United  States 
has  been  developing  an  ASAT  that  can  detect  low  altitude 
satellite  targets  at  a  range  of  a  few  hundred  kilometers. 

If  a  similar  ASAT  is  to  be  effective  against  higher 
orbiting  targets,  the  ASAT  sensor  must  detect  these  targets 
at  a  much  greater  range. 

The  current  United  States  ASAT  uses  an  infrared  sensor 
to  acquire  its  potential  targets  (5).  The  use  of  an 
infrared  sensor  allows  the  sensor  to  acquire  targets 
without  the  need  for  reflected  sunlight  since  any  satellite 
emits  infrared  radiation.  Most  satellites  are  at  a 
temperature  of  275  to  300  degrees  Kelvin.  The  wavelength 
where  the  maximum  energy  output  occurs  for  these 
temperatures  occurs  around  ten  microns. 

Unfortunately,  satellites  are  not  the  only  celestial 
objects  radiating  energy  at  a  wavelength  of  ten  microns. 


Energy  sources  include  planets,  asteroids,  comets,  the 
earth  itself,  interplanetary  and  interstellar  dust  and  gas, 
as  well  as  all  the  stars  in  our  Milky  Way  galaxy.  It  is 
the  vast  number  of  stars  in  the  galaxy  and  the 
interplanetary  dust  within  the  solar  system  that  causes  the 
greatest  source  of  background  radiation  at  ten  microns. 

Detecting  a  target  satellite  against  such  a  background 
is  made  more  difficult  as  the  target  is  moved  further  and 
further  from  the  sensor.  To  complicate  matters  even  more, 
the  celestial  infrared  background  is  by  no  means  uniform. 
The  amount  of  background  radiation  can  vary  by  more  than 
two  orders  of  magnitude.  The  majority  of  stars  are  located 
close  to  the  galactic  equator,  with  fewer  stars  at  higher 
latitudes. 

Two  problems  arise.  What  kind  of  infrared  sensor  is 
required  to  detect  target  satellites  against  such  a 
background?  What  is  the  maximum  distance  the  sensor  can  be 
from  the  target  and  still  expect  to  acquire  its  target? 

This  distance  will  depend  upon  the  complexity  of  the 
infrared  sensor,  the  scan  time,  and  the  area  of  the  sky  the 
target  is  located  in. 

This  thesis  will  determine  values  for  the  infrared 
background  and  will  estimate  the  maximum  range  a  sample 
sensor  can  be  from  a  target,  and  still  detect  it  reliably. 

Chapter  2  of  this  thesis  will  present  background 
material  on  celestial  infrared  radiation  and  infrared 


sensors  used  in  spacecraft. 

Chapter  3  will  discuss  the  two  sources  of  infrared  sky 
data  and  introduce  the  Walker  Sky  Model  to  determine  the 
infrared  radiance  from  different  regions  of  the  sky. 

Chapter  4  will  present  the  parameters  of  sensor  design 
and  determine  the  maximum  detection  ranges  for  a  sample 
sensor . 

Chapter  5  will  present  analysis  of  the  Walker  Sky 
Model,  and  discuss  its  accuracy. 

Chapter  6  will  present  needed  additions  to  the  Walker 
sky  model  and  revised  maximum  sensor  ranges. 

Chapter  7  will  present  conclusions  drawn  from  the 
analysis  and  make  recommendations  for  further  studies. 
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The  backdrop  for  this  problem  lies  in  two  sections 


First,  what  data  exists  on  the  infrared  celestial 
background?  Second,  what  infrared  sensors  have  been  used 
in  the  past  to  detect  exoatmospheric  objects? 

Infrared  Sky  Surveys. 

Two  surveys  of  infrared  celestial  radiation  around  the 
wavelength  of  ten  microns  were  made  during  the  1970's  and 
one  survey  was  made  in  the  current  decade.  The  first  one 
was  the  Celestial  Mapping  Program  (CMP)  by  Grish  and  Lynch 
(23:2-8).  The  second  survey  was  the  Air  Force  Geophysical 
Laboratory  Four  Color  Sky  Survey  by  Price  and  Walker  (27). 
The  Air  Force  Geophysical  Lab  (AFGL)  survey  has  better 
positional  accuracy  than  the  CMP,  but  the  AFGL  survey  has 
gaps  in  its  coverage. 

The  most  recent  survey  of  the  infrared  celestial 
background  was  the  Infrared  Astronomical  Satellite  (IRAS) 
program.  The  IRAS  survey  results  contain  more  than  two 
hundred  thousand  sources  at  wavelengths  greater  than  twelve 
microns.  (2).  Unfortunately,  only  a  small  portion  of  the 
IRAS  data  catalog  is  available.  This  data  is  part  of  the 
IRAS  "Minisurvey"  that  was  conducted  before  the  all-sky 
survey  was  made.  This  minisurvey  consists  of  a  region 
approximately  300  square  degrees  in  size.  This  minisurvey 


was  made  to  test  the  IRAS  components  before  the  main  survey 
began.  The  data  from  this  minisurvey  includes  observations 
of  sources  with  irradiances  values  on  the  order  of  1.0E-18 
W-cm"^-um-1  or  larger  (29). 

Several  classes  of  objects  are  the  most  significant 
contributors  to  the  infrared  sky  background  at  ten  microns. 
These  include  late  type  K  and  M  stars,  stars  deeply 
reddened  by  absorption  from  interstellar  dust,  H-II  regions 
(pronounced  H-two)  (23:2-8),  and  the  zodiacal  light 
particles  (19:458).  H-II  regions  are  interstellar  clouds 
of  ionized  hydrogen  gas.  The  zodiacal  light  is  reflected 
sunlight  from  interplanetary  dust  within  our  solar  system 
(1:253).  Analysis  by  Price  and  Walker  shows  that  over  half 
the  sources  listed  in  the  AFGL  catalog  are  from  either  H-II 
regions  or  stars  reddened  by  interstellar  dust.  The  H-II 
regions  in  the  sun's  neighborhood  have  considerable  angular 
size,  ranging  from  several  arc  minutes  to  several  degrees 
across.  Most  of  the  stellar  and  H-II  sources  are  confined 
within  five  degrees  of  the  galactic  plane.  Walker  and 
Price  also  state  that  many  extended,  low  intensity 
celestial  objects  are  not  included  in  the  AFGL  catalog 
because  of  their  low  surface  radiance  (23:2-8). 

The  measurement  of  emission  from  the  interplanetary 
dust  that  causes  the  zodiacal  light  has  been  more  difficult 
because  of  the  large  angular  extent  of  the  emission  source. 
Briotta,  Piper,  and  Houck  made  observations  in  the  eight  to 


twelve  micron  range,  but  their  observations  were  seriously 
contaminated  with  radiation  from  the  earth  (24:771). 
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Observations  have  also  been  made  by  the  AFGL,  but  the  data 
is  measured  relative  to  the  non-ecliptic  background,  and  is 
not  in  terms  of  absolute  radiances  (23:2-10). 

In  addition  to  the  three  classes  of  ten  micron 
radiation  mentioned  above,  asteroids  also  present  a 
possible  background  source.  The  total  number  of  asteroids 
in  the  solar  system  is  unknown,  but  is  estimated  to  be  far 
in  excess  of  the  3000  for  which  orbits  have  been 
calculated.  From  observations  at  eleven  microns,  asteroids 
temperatures  have  been  estimated  to  be  Detween  200  to  260 
degrees  Kelvin.  This  is  the  range  in  which  earth  orbiting 
satellites  would  be  found  (23:2-11).  Murdock  (21)  made  a 
study  to  predict  the  number  of  asteroids  visible  an  any  one 
time.  The  asteroids  are  concentrated  within  seven  degrees 
of  the  ecliptic  plane.  The  density  of  asteroids  along  the 
ecliptic  plane  may  be  several  per  square  degree  and 
contribute  1.0E-19  W-cnT^-unT1  to  the  background  irradiance 
at  eleven  microns  (23:2-11). 

A  recent  work  by  Grasdalen  et  al  (9:413)  tried  to 
correlate  AFGL  sky  survey  sources  to  ground  based 
observations.  The  effort  was  to  isolate  that  portion  of 
the  AFGL  sky  survey  that  included  only  stellar  objects. 
Another  objective  of  this  study  was  to  estimate  the  type 
and  number  of  sources  the  Infrared  Astronomical  Satellite 
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(IRAS)  would  be  likely  to  encounter.  The  study  estimated 
that  IRAS  would  see  about  100,000  more  sources  than  the 
AFGL  survey,  but  the  mixture  of  sources  would  not  be 
substantially  different  than  those  seen  by  the  AFGL  survey. 

As  mentioned  previously,  the  most  recent  work  in  the 
measurement  of  the  celestial  background  between  eight  and 
fifteen  microns  was  conducted  by  IRAS.  Preliminary 
analysis  of  data  indicates  an  additional  ring  of  dust 
particles  further  from  the  sun  than  the  ring  that  causes 
the  zodiacal  light.  The  temperature  of  the  dust  is 
estimated  to  be  between  165  and  200  degrees  Kelvin  (28:5). 
Some  of  data  gathered  from  earlier  investigations  was 
confirmed  by  IRAS.  This  included  the  contribution  to  the 
background  from  very  large  emission  sources  (H-II  regions) 
located  within  the  galactic  plane,  and  the  fact  that  the 
galactic  nucleus  was  shown  to  be  an  especially  strong 
emitter  (2:122).  However,  the  IRAS  measured  zodiacal 
emission  measurements  exceeds  previously  published 
measurements  at  eleven  microns  (11:617).  The 
interpretation  of  the  data  indicates  that  the  emission 
stems  from  optically  thin  clouds  of  dust  which  are  not  at 
a  single  temperature. 

Infrared  SpasLe.  Sensors 

Various  infrared  space  sensors  have  been  designed  or 
developed  for  detecting  satellites  over  the  past  decade. 
These  can  be  divided  into  two  relevant  categories:  ground 
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based  and  space  based. 

Work  in  this  area  for  a  deep  space  tracking  network 
was  accomplished  in  the  early  1970's.  This  work  was  done 
before  the  decision  was  made  to  use  a  visual  sensor  for  the 
Ground-Based  Deep  Space  Surveillance  System  (GEODSS). 
Proposals  included  an  airborne  Long  Wave  Infrared  (LWIR) 
sensor  deployed  on  aircraft  above  thirty  thousand  feet 
(34).  NASA's  Kuiper  airborne  infrared  observatory  has  a  C- 
141  and  uses  a  similar  approach  for  astronomical  research. 
The  most  successful  ground  based  development  in  this  area 
is  the  Maui  Optical  Tracking  and  Identification  Facility 
(MOTIF)  on  Mount  Haleakala  in  Hawaii.  This  sensor 
continues  to  do  active  research  and  spacetrack 
augmentation. 

The  MOTIF  sensor  consists  of  a  1.2  meter  Cassegrain 
telescope  with  an  LWIR  radiometer  attached.  The  sensor 
array  consists  of  25  germanium-cadmium  detectors  cooled  to 
13  degrees  Kelvin.  The  infrared  range  is  from  3  to  22 
microns  in  seven  different  bands.  The  field  of  view  for 
the  sensor  is  80  arc  seconds. 

The  MOTIF  sensor  has  been  used  to  determine  LWIR 
signatures  of  foreign  satellites  as  well  as  investigations 
into  the  LWIR  radiance  of  the  atmosphere.  The  objectives 
of  the  satellite  infrared  signature  program  are  to  obtain 
both  passive  and  active  infrared  signatures  in  the  8  to  13 
and  17  to  21  micron  bands.  This  program  tries  to  gather 


sufficient  data  on  specific  satellites  to  allow  the 
prediction  of  signatures  for  satellites  that  cannot  or  have 
not  been  observed  (30). 

In  the  1974-5  timeframe,  several  studies  were 
commissioned  to  design  a  space  based  satellite  tracking 
satellite.  Studies  were  made  by  Philco  Ford  Corporation, 
Lockheed  Missiles  and  Space  Corporation,  and  Rockwell 
International  (15). 

After  the  early  studies  were  completed,  an 
experimental  satellite  was  contracted  for  in  order  to 
increase  the  data  base  necessary  to  accurately  design  an 
LWIR  space  surveillance  system  satellite.  This 
experimental  satellite  was  known  as  the  Space  Infrared 
Experiment  (SIRE).  The  satellite  was  originally  scheduled 
to  be  launched  in  1981,  but  the  project  was  delayed  and 
eventually  abandoned  because  of  large  anticipated  costs 
(33).  However,  the  technology  developed  for  this  program 
may  be  used  on  a  future  space  surveillance  system  (35). 

The  infrared  sources  the  satellite  was  to  investigate 
included  earth  limb  radiance,  zodiacal  light,  stellar 
radiation,  and  asteroids,  as  well  as  satellites.  The 
satellite  was  designed  to  determine  which  spectral  band  in 
the  infrared  gave  the  best  possibility  of  detecting  a 
target  satellite.  Also,  this  satellite  would  have  provided 
a  test  of  the  optics  and  two  dimensional  detector  array 


The  SIRE  satellite  was  to  be  a  cryogen  cooled  sensor 
covering  the  7  to  23  micron  range  in  five  separate  bands. 

It  was  to  use  a  three  mirror  anastigmat  optical 
configuration  with  a  collecting  area  of  435  square 
centimeters.  This  arrangement  would  allow  for  detection 
and  tracking  of  objects  in  both  low  earth,  and 
geosynchronous  orbits  (18). 

Several  spacecraft  have  used  long  wavelength  infrared 
sensors  as  part  of  their  instrumentation.  These  spacecraft 
include  Landsat  D,  the  Defense  Meteorological  Satellite 
(DMSP),  the  United  States  ASAT  test  project,  and  the 
Infrared  Astronomical  Satellite  (IRAS)  (11). 

The  IRAS  sensor  was  used  to  conduct  a  survey  of  the 
entire  sky  from  an  altitude  of  900  kilometers.  This  sensor 
conducted  observations  in  four  infrared  bands:  12,  25,  60, 
and  100  microns.  The  detectors  for  these  bands  were 
silicon-arsenide,  silicon-tin,  germanium-gallium,  and 
germanium-gallium  respectively.  The  entire  sensor 
consisted  of  a  telescope  0.57  meters  in  diameter  with  a 
five  meter  focal  length.  The  focal  plane  array  with  the 
detectors  was  cooled  to  less  than  three  degrees  Kelvin. 

The  two  mirror  beryllium  optical  system  (Ritchey-Chretien) 
was  cooled  to  less  than  ten  degrees  Kelvin  by  liquid 
helium. 

The  focal  plane  assembly  was  constructed  so  that  every 
infrared  source  crossing  the  60  arc  minute  field  of  view 


could  be  seen  by  at  least  two  detectors  in  each  wavelength 
band.  The  position  of  each  source  seen  by  the  16  detectors 
in  the  12  micron  band  lies  in  an  uncertainty  ellipse 
approximately  45  x  9  arc  seconds.  Each  detector  was 
covered  by  a  field  lens  and  a  filter  giving  a  mean  in-band 
transmission  of  83  percent  for  the  12  micron  band.  The 
system  in-band  noise  equivalent  flux  density  (NEFD)  was  on 
the  order  of  1.0E-19  W-cm-^  averaged  across  the  spectral 
bandpass.  (20:69). 
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III.  The 


There  are  two  major  catalogs  of  infrared  stars.  One 
is  from  the  Air  Force  Geophysics  Lab  Sky  Survey.  This  data 
catalog  was  compiled  from  nine  rocket  flights  between  1971 
and  1974.  The  catalog  was  published  in  1976.  At  eleven 
microns,  90  percent  of  the  sky  was  surveyed  at  least  once 
and  about  60  per  cent  of  the  sky  was  mapped  two  or  more 
times  (27:9).  The  AFGL  catalog  lists  1151  sources  at 
eleven  microns.  Through  data  reduction  and  spurious  source 
elimination,  the  AFGL  considers  the  catalog  complete  down 
to  an  irradiance  level  of  2.2E-16  W-cm~2-um-1  at  eleven 

microns.  The  catalog  does  list  some  sources  as  faint  as 
1.2E-16  W-cm-2-um-1  (27:14).  These  1151  sources  are  not 
spread  out  evenly  across  the  sky,  but  the  majority  are 
within  plus  or  minus  20  degrees  of  the  galactic  plane. 

The  second  catalog  is  the  NASA  Catalog  of  Infrared 
Observations.  This  catalog  is  a  listing  of  all  the 
astronomical  observations  in  the  infrared  between  1  and 
1000  microns,  which  were  published  in  the  scientific 
literature  from  1965  through  1982.  This  catalog  contains 
over  85,000  individual  observations  of  about  10,000 
different  infrared  sources.  The  catalog  was  published  in 
May  1984  (8:1).  This  catalog  contains  the  data  from  the 
AFGL  catalog,  and  additional  data.  The  NASA  catalog  has 
objects  listed  as  faint  as  1.2E-18  W-cm_2-um“1,  but  these 


entries  are  sporadic  and  far  from  complete. 

The  two  major  infrared  source  catalogs  only  contain 
reliably  complete  data  as  faint  as  1.0E-16  W-cm-2-um“*. 

Since  spacecraft  target  intensities  can  be  as  faint  as 
1.0E-18  W-cm-2-um-l  (10:61),  a  source  of  data  for  stars  at 
this  irradiance  level  and  lower  must  be  found.  IRAS  data 
includes  stars  this  faint,  but  its  entire  data  catalog  is 
unavailable  at  present. 

The.  Walker  Sky  Model 

Walker  has  developed  a  model  for  estimating  the 
density  and  hence  irradiance  levels  for  infrared  sources 
within  the  galaxy  (36:3-28).  This  model  estimates  the 
number  of  sources  observed  at  irradiance  levels  from  1.0E- 
14  to  1.0E-20  W-cm-2-um_1  at  various  galactic  latitudes 

and  longitudes. 

The  model  represents  the  galaxy  as  a  planar  disk  with 
a  dense  central  bulge.  The  center  of  the  galaxy  is  at  a 
distance  R  =  1.0E4  parsecs  (1  parsec  =  3.09E18  cm)  from  the 
Sun.  Galactic  coordinates  are  1  for  longitude,  measured 
from  the  Sun  galactic  center  line,  and  b  for  latitude, 
measured  from  the  galactic  disk.  The  galaxy  is  shown 
graphically  in  figure  1. 

The  model  of  infrared  sources  is  considered  to  be  a 
Gaussian  distribution  of  the  form  A  exp  -  (X  /(7)2  where  X 
is  the  distance  above  the  galactic  plane  and  (J  is  the 


*  (arbitrary  point) 


standard  deviation  of  the  galactic  distribution  of 
infrared  sources  in  terms  of  perpendicular  distance  from 
the  galactic  plane  (10).  Several  classes  of  objects 
contribute  to  the  infrared  source  model.  These  are  dwarf 
stars  (mass  less  than  10  solar  masses,  1  solar  mass  = 
1.991E33  grams),  giant  stars  (mass  less  than  50  solar 
masses),  supergiant  stars  (mass  greater  than  50  solar 
masses),  and  HII  regions  (1).  The  dispersions,  (J^  ,  normal 
to  the  galactic  plane  are  160  parsecs  for  dwarf  stars,  450 
parsecs  for  giant  stars,  100  parsecs  for  supergiant  stars, 
and  100  parsecs  for  HII  regions. 

The  nucleus  of  the  galaxy  contains  more  sources  per 
volume  than  does  the  rest  of  the  galaxy.  The  Gaussian 
distribution  is  modified  by  the  factor. 


where  Rc  is  the  radius  of  the  galactic  nucleus  =  9.27E21 
cm,  and  R^  is  the  distance  (shown  in  figure  1)  in  the 
galactic  plane  from  the  galactic  center  to  some  arbitrary 
point. 

This  model  is  represented  by  a  single,  though 
complex  equation: 


M 


where 


N  =  the  number  of  infrared  sources  per  square  degree  in 
some  arbitrary  direction  along  by  r^  (specified 
by  angles  b,l),  producing  irradiance  at  the 
observer  which  is  greater  than  E^  , 

E  =  the  limiting  spectral  irradiance  in  W-cm"^-um-^, 

A  ' ’**  c 

I1  =  a,  X1  #  the  mean  spectral  radiant  intensity  of 
*■ 

the  ith  class  of  celestial  objects  in  watts/ster, 
n^  =  the  mean  space  density  of  the  ith  class  of 
sources  in  numbers/cm^, 

1  =  galactic  longitude  in  degrees  east  from  a  line 

between  the  observer  and  galactic  center, 
b  =  galactic  latitude  in  degrees  above  or  below  the 
plane  of  the  galaxy  as  seen  from  the  observer's 
location, 

R0  =  the  distance  from  the  Sun  (and  observer)  to  the 
galactic  center  =  3.09E22  cm, 

Rc  =  the  radius  of  the  galactic  nucleus  =  9.27E21  cm, 

-  the  dispersion  of  the  ith  class  normal  to  the 
galactic  plane  in  cm, 

X  =  the  wavelength, 

R*  =  I.  .  e"'  cos(b)?+  R*  -2R  i'*  E^cos(b)  cos(l) 

+  I.  E  cos  (b) , 

aifCi  =  constants  observationally  determined  which 

describe  the  intensity  of  the  generic  member  of 
the  ith  class,  and 
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k 


=  1.52E-4  ster-deg-^. 

The  values  for  (Jk  ,  n^,  a^,  and  for  the  model  are 
derived  from  sky  survey  and  other  astronomical  data 
(10:121)  and  are  given  in  table  I.  The  irradiances  caused 
by  emission  nebulae  (H-II)  have  been  included  with  the 
supergiants  because  they  have  similar  space  distributions. 

The  value  for  is  determined  using  the  law  of 
cosines  (see  figure  1) 

R*  =  R*  +  (r  cos(b)  )J  -  2R<J(rtcos(b))  cos(l), 
and  the  point  source  intensity-irradiation  relation 


The  equation  for  the  model,  while  not  obvious,  is  a 
best  fit  to  the  AFGL  and  other  data  (10),  and  does  check 
out  to  be  dimensionally  correct.  A  discussion  of  its 
construction  is  found  in  Appendix  D. 


Table  I 

Constants  for  the  Galactic  Source  Model 


Class 

Vi 

ni 

ai 

ci 

Dwarfs 

4.94E+20 

2.63E-56 

2.77E+24 

-3.693 

Giants 

1.32E+21 

9.46E-60 

1.07E+28 

-1.915 

H  II  regions 

3.09E+20 

8.07E-62 

1.09E+29 

-1.108 

&  supergiants 

units 

cm 

-3 

stars-cm 

-)  -  » 
W-sr  -urn 

— 
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The  model  will  yield  the  number  of  sources  for  any 
particular  direction  defined  by  b  and  1,  and  for  a 
particular  wavelength  defined  by  X .  To  evaluate  the 
integral  for  each  wavelength  region,  the  interval  of 
integration  was  broken  into  three  latitude  ranges  and 
numerical  integration  was  done  by  means  of  16-point 
Legendre-Gauss  quadrature.  This  was  done  by  Walker  (36:3- 
30)  and  Greenebaum  (10:127). 

The  results  the  integration  for  eleven  microns  are 
shown  in  table  II  (10:134).  Since  the  galaxy  is  considered 
to  be  symmetric,  the  values  for  negative  latitudes  will  be 
the  same  as  those  for  positive  latitudes.  The  same  can  be 
said  for  longitude  between  0  and  180  degrees,  and  those 
between  180  and  360  degrees. 

This  table  lists  the  number  of  infrared  sources  per 
square  degree  at  a  wavelength  of  eleven  microns  for 
selected  areas  of  the  galaxy.  The  spectral  irradiance  is 
in  Watts-cm-2-um-1.  Each  item  in  the  body  of  the  table  is 
the  sum  of  all  sources  to  the  limit  of  the  spectral 
irradiance  for  that  row  in  the  table.  For  example,  for 
galactic  latitude  0  degrees,  longitude  0  degrees  there  are 
9400  infrared  sources  at  a  spectral  irradiance  level  of 
1.0E-18  Watts-cm”2-um”^  and  brighter. 

To  find  the  total  radiance  per  square  degree  for  any 
particular  area,  the  contribution  of  each  spectral 
irradiance  level  must  be  calculated.  Since  each  entry  in 
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Irradiance  is  in  Watts-cm 


A.. 

1 


the  table  represents  a  sum  of  all  the  sources  to  that  level 
of  brightness/  the  sources  brighter  than  that  particular 
level  must  be  subtracted  in  order  to  obtain  the  irradiance 
from  that  level  only.  That  isf  to  find  the  number  of 
sources  only  at  an  irradiance  range  of  1.0E-15  Watts-cm  - 
um-1,  to  10E-15  the  number  of  sources  at  1.0E-14  must  be 
subtracted  from  the  number  for  1.0E-15  listed  in  the  table. 
The  number  of  sources  at  a  particular  irradiance  range 
times  the  average  irradiance  of  that  range  will  give  that 
ranges's  contribution  to  the  total  radiance  from  a  galactic 
region.  The  sum  of  the  irradiance  level  contributions 
gives  the  total  radiance  for  that  area. 

For  example,  to  find  the  total  radiance  from  a  region 
around  galactic  latitude  10  degrees,  longitude  45  degrees, 
the  total  number  of  individual  components  must  first  be 
calculated,  as  shown  below. 


From  table  II, 

Irradiance  E^ 
W-cm  -urn 

1.0E-14 

1.0E-15 

1.0E-16 

1.0E-17 

1.0E-18 

1.0E-19 


1.0E-20 


Number  of  Sources 

With  Ea  >  E>e 

4.9E-4 

7.3E-3 

9.7E-2 


2.8E0 

4.7E1 

7.7E1 

7.9E1 


Number  of  Sources  in 


Irradiance  Ranges  of  E 

A o 

Irradiance  Range 


>  1.0E- 

■14 

4.9E-4 

-  0  = 

4.90E-4 

1.0E-14 

to 

1.0E-15 

7.3E-3 

-  4.9E-4 

= 

6.81E-3 

1.0E-15 

to 

1.0E-16 

9.7E-2 

-  7.3E-3 

= 

8.97E-2 

1.0E-16 

to 

1.0E-17 

2.8E0 

-  9.7E-2 

= 

2.70E0 

1.0E-17 

to 

1.0E-18 

4.7E1 

-  2.8E0 

= 

4.42E1 

1.0E-18 

to 

1.0E-19 

7.7E1 

-  4.7E1 

= 

3.00E1 

1.0E-19 

to 

1.0E-20 

7.9E1 

-  7.7E1 

— 

2.00E0 

Next,  the  individual  numbers  must  be  multiplied  by  the 
irradiance  levels  to  find  the  contribution  from  that 
irradiance  level.  That  is: 


Number 

Irradiance 
W-cm“2-um  1 

(4.90E-4) 

X 

(1.0E-14)  = 

(6.81E-3) 

X 

(1.0E-15)  = 

(8.97E-2) 

X 

(1.0E-16)  = 

(2.70E0) 

X 

(1.0E-17)  = 

(4.42E1) 

X 

(1.0E-18)  = 

( 3.00E1) 

X 

(1.0E-19)  = 

(2.00E0) 

X 

(1.0E-20)  = 

The  irradiance  at  the  botto 


Contribution 

4.90E-18 
6.81E-18 
8.97E-18 
2.70E-17 
4.42E-17 
3.00E-1 8 
2.00E-20 

n  of  each  range  is  chosen 


as  the  average  value  for  the  range  because  the  number  of 
sources  is  greatest  at  lower  irradiances. 


Since  table  II  gives  the  number  of  sources  per  square 
degree  in  each  irradiance  range,  summing  the  radiance 
contribution  from  each  irradiance  range  will  give  the  total 
radiance  expected  from  a  particular  1  degree  by  1  degree 
region  of  the  sky. 

Summing  the  individual  contributions  calculated  above 
gives  the  total  radiance  for  the  1  degree  by  1  degree 
region  of  the  sky  around  latitude  10  degrees,  longitude  45 
degrees.  It  is  9.49E-17  Watts-cm-2-um“l-deg“2,  or  6.245E- 
13  Watts-cm_2-um~*-ster“^. 

The  radiance  calculations  for  each  sky  area  of  table 
II  are  shown  in  table  III.  Table  III  shows  total  radiance 
in  Watts-cm-2-um“1-ster“1  for  the  galactic  positions  shown. 
The  radiances  calculated  in  this  table  are  in  general 
agreement  with  those  found  in  the  Infrared  Handbook  (36:3- 
32).  From  this  table  a  contour  map  can  be  created  showing 
the  expected  radiances  as  a  function  of  galactic  position. 
This  is  shown  in  figure  2. 
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TOTAL  RADIANCE  IN  WATTS-CM_2-UM-1-STER“1 


GALACTIC 


GALACTIC  LONGITUDE 


■TITUDE 

0 

45 

90 

180 

0 

3.31E-11 

1.00E-11 

4.81E-12 

2.97E-12 

2 

1.25E-11 

3.76E-12 

1.95E-12 

1.28E-12 

5 

2.43E-12 

1.40E-12 

9.52E-13 

7.18E-13 

10 

7.19E-13 

6.25E-13 

5.07E-13 

4.20E-13 

20 

2.88E-13 

2.75E-13 

2.51E-13 

2.28E-13 

30 

1.91E-13 

1.82E-13 

1.73E-13 

1.62E-13 

45 

1.24E-13 

1.22E-13 

1.18E-13 

1.12E-13 

60 

9.61E-14 

9.55E-14 

9.36E-14 

9.11E-14 

90 

8.08E-14 

8.08E-14 

8.08E-14 

8.08E-14 

Lgure  2.  Radiance 


IV.  SENSOR  0££I£H 


Various  combinations  of  infrared  detectors  can  be 
used  to  make  up  a  sensor.  Sensors  can  use  a  single 
detector,  a  linear  array  of  detectors,  or  a  two  dimensional 
array  of  detectors  as  in  a  charge-coupled  device  (CCD). 

The  matrix  of  detector  cells  can  be  spun  around  an  axis  to 
create  a  scanning  sensor  or  may  be  stationary  to  create  a 
staring  array.  The  detector  matrix  can  be  held  stationary 
and  a  reticle  moved  in  front  of  the  detectors  to  cause  some 
angular  or  positional  modulation  of  the  incoming  infrared 
radiation. 

Detector  cells  for  infrared  light  at  around  ten 
microns  can  be  made  of  several  combinations  of  materials, 
the  most  common  being  mercury-cadmium-telluride  or  silicon- 
arsenide.  Silicon-arsenide  was  used  on  the  IRAS  sensor 
(20).  Mercury-cadmium-telluride  is  used  in  the  LANDSAT  D 
thematic  mapper  and  DMSP  (12).  Detector  sensitivity  varies 
with  the  material,  but  all  must  be  cooled  to  cryogenic 
temperatures  to  be  useful  in  a  infrared  sensor.  This 
sensitivity,  the  quantum  efficiency,  is  the  average  number 
of  electrons  produced  per  photon  hitting  the  detector 
material.  The  quantum  efficiency  varies  from  about  10  to 
50  per  cent  (19:1) . 

In  order  for  an  ASAT  to  successfully  attack  a  target 
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spacecraft,  the  ASAT's  sensor  must  acquire  the  target  at  a 
large  enough  distance  so  that  orbital  corrections  can  be 
made  for  intercept.  The  sensor's  design  must  be  sensitive 
enough  to  ensure  reliable  target  detection,  but  not  so 
complex  that  it  becomes  too  expensive  for  a  single  use  type 
mission. 

Besides  detector  sensitivity,  the  spatial  resolution 
of  the  sensor  is  also  important.  This  is  the  ability  to 
discriminate  between  two  closely  space  objects.  The 
resolution  is  a  function  of  the  size  of  the  collecting 
optics  the  sensor,  the  detector  size  and  spacing  in  the 
array,  atmospheric  seeing,  and  the  wavelength  of 
operation.  For  a  sensor  whose  resolution  is  diffraction 
limited  by  its  optics,  the  resolution  Q 


=  1.22 


(4.1) 


where 

X  =  the  operating  wavelength  and 
D  =  the  diameter  of  the  optics, 
is  the  angular  distance  in  radians  between  the  closest 
two  objects  that  will  be  discernible  as  two  objects  instead 
of  one  object. 

Resolution  is  not  solely  dependent  upon  optics, 


however;  the  number  of  detectors  and  their  size  determines 
the  spatial  resolution  as  well.  If  two  point  sources  are 


to  be  detected  as  two  point  sources,  then  their  images  must 
fall  on  two  different  detectors.  The  minimum  resolvable 
angle  Q 

d 

Q  - -  (4.2) 

F 

where 

d  =  the  spacing  between  detectors, 

F  =  the  focal  length  of  the  optic  system, 

0  =  the  resolvable  angle,  in  radians. 

As  an  example,  the  infrared  detector  on  LANDSAT  D  has 
a  resolution  of  about  40  meters  from  a  distance  of  800 
kilometers.  If  this  same  resolution  existed  for  an  ASAT 
sensor,  the  ASAT  could  distinguish  between  two  targets  a 
half  a  kilometer  apart  at  a  range  of  10,000  kilometers. 

The  ASAT  should  not  need  a  resolution  capability  this  fine 
however,  because  geosynchronous  spacecraft  are  kept  at 
least  600  kilometers  apart  (17). 

Since  the  ASAT's  sensor  is  expected  to  acquire  the 
target  at  several  hundred  to  several  thousand  kilometers 
distance,  any  target  will  be  only  a  point  source  of 
infrared  radiation.  The  assumed  sensor  will  be  operating 
with  a  wavelength  bandpass  around  ten  microns  where  most 
spacecraft  with  a  temperature  of  250  to  350  degrees  K. 
radiate  the  most  energy.  The  assumed  target  spacecraft 


will  radiate  energy  with  intensity  that  is  function  of  the 
emissivity  of  the  surface  and  area  facing  the  sensor. 


In  general,  real  spacecraft  are  not  isotropic 
radiators.  The  amount  of  energy  radiated  is  dependent  upon 
the  viewing  angle  and  all  the  emissivities  of  the 
spacecraft  outer  components.  Each  different  type  of 
surface  material  -  solar  cells,  antennas,  insulation, 
paint,  etc.  -  has  its  own  emissivity.  It  is  very  difficult 
to  predict  the  radiant  characteristics  of  a  particular 
target  spacecraft  whose  projected  area  is  constantly 
changing  with  respect  to  the  viewing  angle.  Infrared 
irradiances  for  generic  satellites  at  a  temperature  T  are 
best  described  in  terms  of  equivalent  emissivity-areas. 

In  order  to  design  a  sensor  for  maximum  efficiency, 
several  factors  must  be  considered: 

the  collecting  area  of  the  sensor, 
the  quantum  efficiency,  that  is,  the  number  of 
measurable  electrons  generated  per  photon  input  to 
detector, 

the  percentage  of  radiation  transmitted  through 
the  sensor  optics, 

the  spectral  bandpass  of  the  detector, 
the  angle  of  incidence  of  the  radiation  onto  the 
collecting  area  of  the  detector, 
the  field  of  view  of  the  detector  in 
steradians,  and 
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Ts  -  the  interval  between  samples  of  the  detector 
output  current. 

The  design  of  the  sensor  will  determine  all  of  these 
factors.  One  of  the  simplest  sensors  is  a  single  detector 
with  a  spinning  reticle.  This  is  placed  at  the  focal  plane 
of  the  sensors  optical  system.  The  reticle  consists  of  a 
circular  mask  with  a  single  rectangular  slit.  A  possible 
layout  is  shown  in  figures  3  and  4.  The  size  of  the 
reticle  slit  will  determine  how  much  of  the  detector  is 
illuminated  at  any  one  time.  When  the  reticle  is  spun 
about  its  center,  objects  appearing  near  the  center  of  the 
reticle  remain  visible  for  a  longer  period  of  time  than  do 
objects  near  the  edge  of  the  reticle.  It  is  this  period  of 
visibility  coupled  with  the  time  that  the  object  is  seen  in 
the  rotation  period  that  will  determine  the  direction  of 
the  object  relative  to  the  axis  of  the  sensor. 

The  output  of  the  detector  must  be  sampled  often 
enough  so  the  ASAT  guidance  system  can  make  corrections  in 
the  trajectory  if  needed.  The  speed  of  the  reticle 
rotation  must  be  synchronized  with  the  sampling  rate  so 
there  are  no  gaps  in  coverage  of  the  detector,  yet 
positional  accuracy  is  the  best  possible. 

Since  the  size  of  the  reticle  slit  determines  how  much 
of  the  total  field  of  view  the  detector  will  see  at  any 
particular  instant,  it  must  be  taken  into  account  when 


designing  a  possible  ASAT  sensor.  That  is,  if  the  area  of 
the  reticle  slit  is  l/10th  of  the  total  area  of  the 
detector,  the  resulting  instantaneous  field  of  view  of  the 
sensor  is  only  l/10th  of  the  total  field  of  view  of  the 
optics . 

The  sampling  time  determines  the  positional  accuracy 
of  the  sensor.  The  more  often  the  output  current  is 
sampled,  the  finer  resolution  the  sensor  has.  For  example, 
a  sensor  with  a  reticle  spinning  ten  rotations  per  second 
and  reticle  slit  20  degrees  wide,  the  field  of  view  is 
divided  into  18  different  sections.  If  the  current  is 
sampled  four  times  during  the  time  it  takes  the  reticle  to 
move  twenty  degrees,  the  current  is  sampled  once  every  Ts  = 
1.39E-3  seconds.  The  sampling  rate  is  the  inverse  of  this, 
or  about  720  times  a  second. 

Since  the  output  current  is  sampled  every  5  degrees  of 
reticle  rotation,  the  positional  accuracy  relative  to  the 
circular  scan,  is  1/72  of  the  total  scan.  This  breaks  the 
field  of  view  into  72  circular  sectors.  The  positional 
accuracy  relative  to  the  distance  along  a  radius  of  the 
detector  face  is  more  difficult  to  determine.  Since  the 
detector  slit  is  rectangular,  the  part  of  the  detector 
closest  to  the  center  is  illuminated  longer  than  portions 
near  the  outer  edge.  This  means,  if  a  target  is  close  to 
the  center  of  the  optical  axis,  it  will  remain  visible  for 
more  current  samples  than  if  the  target  was  near  the  edge 


of  the  field  of  view.  For  this  sensor,  targets  at  the  edge 
of  the  field  of  view  will  be  visible  for  4  sampling  periods 
while  objects  at  the  center  will  be  visible  for  32  sampling 
periods.  This  will  divide  the  field  of  view  into  28 
different  radial  zones.  The  radial  zones  have  more  area 
toward  the  edge  of  the  detector  and  less  area  toward  the 
center.  The  circular  sectors  all  have  the  same  area. 

It  is  the  intersection  of  these  sectors  and  zones 
that  form  the  resolution  element  (see  figure  5).  The  size 
of  any  resolution  element  is  determined  by  its  distance 
from  the  center  of  the  detector.  At  a  point  half  way  out 
from  the  center,  a  target  would  be  sampled  8  times  per 
reticle  revolution  (see  figure  5).  The  size  of  the  radial 
zone  having  8  samples  per  reticle  field  of  view  is  about 
17%  of  the  total  field  of  view.  Since  the  field  of  view  is 
also  broken  up  into  72  circular  sectors,  multiplying  0.17 
times  1/72  will  give  the  part  of  the  total  field  of  view 
contained  in  this  particular  resolution  element.  This  is 
about  0.0024.  Since  the  total  field  of  view  is  about  78.5 
square  degrees,  this  resolution  element  is  about  0.434 
degrees  on  a  side. 

The  object  of  the  ASAT  sensor  design  is  to  determine 
how  far  away  a  proposed  sensor  can  be  from  its  target  and 
still  detect  it  reliably.  The  target  must  be  detected 
against  the  infrared  celestial  background  discussed  in  the 
previous  section.  To  do  this  a  series  of  calculations  must 


be  performed.  To  perform  these  calculations,  a  set  of 
parameters  for  the  proposed  sensor  must  be  used.  They 
involve  the  signal-to-noise  ratio,  the  signal  current 
produced  in  the  sensor  from  a  target,  and  the  signal 
current  produced  in  the  sensor  from  the  celestial 
background.  Table  IV  lists  a  set  of  possible  sensor 
parameters . 

TABLE  IV 

Example  ASAT  Sensor  Parameters 

Optics  Diameter  *  .2  meters,  S  =  D2  =  0.0314  m2 

4 

Optics  Cooled  to  10°  Kelvin 
Optics  Transmission  -  =  0.9 

O 

Field  of  View  (Half  Angle)  =  5  degrees 

Solid  Angle  0  -  7T  6  -  0.024  Steradians 

Fov 

Detector:  Si:As 

NEFD  =  1.0E-18  Watts-cm-2 

Focal  Plane  Cooled  to  3°  Kelvin 

Cryogen:  Liquid  Helium 

Detector  Quantum  Efficiency  -  TJ(xj  a  .25 

Detector  Sample  Rate  =  720  samples  per  second 

Reticle  Spin  Rate  =  10  rotations  per  second 

Reticle  Slit  Area  =  11%  of  total  detector 

Wavelength  Window  =  9  to  13  microns 

Average  Wavelength  X  =  11  microns 

Desired  Signal-to-Noise  Ratio  -  (S/N)  =  5 


From  a  point  source  target  the  amount  of  current  iT 
output  from  a  detector  (12)/ 


e  Sft  cos  6a. 
her1 


At 


^2 


(4.3) 


where 

1^  =  Target  intensity  in  Watts/ster 
e  =  charge  on  the  electron/  1.6E-19  coulombs 
h  =  Planck's  constant,  6.62E-34  joule-sec 
c  =  speed  of  light,  3.0E8  m/sec 
r  =  range,  the  distance  from  the  target  to  the 
sensor 

From  an  extended  background  the  amount  of  current  iB 
output  from  a  detector  is  (12) 

i&-  I  L.M'nmji) a<**  u-4> 

A, 

where  is  the  solid  angle  of  the  filed  of  view 

and  L  is  the  background  radiance  in  Watts-cm“^-um-1-ster  . 
To  achieve  reliable  detection  of  a  target  the  signal- 


to-noise  ratio  (S/N)  must  be  high  enough  so  the 
intermittent  noise  and  background  inhomogeneities  do  not 
frequently  produce  false  targets.  A  signal-to-noise  ratio 


of  five  can  guarantee  few  false  alarms.  The  signal-to- 
noise  ratio  is  defined  as  the  average  value  of  the  current 
just  from  the  target  divided  by  the  variance  of  the  total 
current  from  the  detector.  The  total  current  includes  the 
target  current,  background  current,  and  noise  current. 
Specifically  for  a  sensor  using  current  for  detection, 
where  the  only  source  of  noise  is  photon  noise,  the  signal- 
to-noise  ratio  is  (37:541), 


_ £t _ 

((ir  *  4  i-o )  3  e  Af  )  ^ 


(4.5) 


where 

iT  =  total  current 

iB  =  background  current 

iD  =  dark  or  noise  current 

Af  *  frequency  bandwidth  of  the  circuitry 

e  =  charge  of  the  electron 

To  design  the  sensor  the  analysis  will  proceed  as 
follows: 

1.  Solve  the  Signal-to-Noise  equation  (eq.  4.5)  for 
the  Target  Current  iT. 

2.  Solve  the  Target  Current  equation  (eq.  4.3)  for 


the  Range. 

3.  Determine  the  Background  Current  from  the  Background 


Current  equation  (eq.  4.4). 


j;iia 

% 


4.  Determine  the  value  of  iT  obtained  from  the  Signal- 
to-Noise  equation  (eq.  4.5). 

5.  Substitute  the  value  of  iBto  the  Range  equation 
(eq.  4.3). 

6.  Solve  for  the  range  r. 

1 

Since  f  is  approximately  equal  to  -  ,  where 

2  ts 

ts  the  sampling  time  of  the  detector,  the  signal-to-noise 
ratio  (12),  becomes: 


Solving  for  iT  yields: 
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The  target  current  from  an  irradiance  X  (\)  is: 

A 


Ai 


Evaluating  this  across  a  small  wavelength  window  near 
the  peak  of  the  blackbody  curve  gives: 

i  .  g-Sa.'w).  i  a)  t  a)  7/a)aa 
T '  hc/\ 

This  assumes  Y  (X)  /  7^(7)  are  constant  across 

the  window.  Since  the  sensor  is  outside  of  the  atmosphere 
fAa)  - 1-°-  On  IRAS  7a(\)  was  about  0.9,  which  is  the 
value  assumed  here  (20).  This  leads  to: 


Lj  - 


tli.  1  ^ _  i  (i)  iA(x)  ax 


Kc,  rx 


The  irradiance  from  a  point  source  is: 

Et:  7iJ  --  7‘ 


A, 


Substituting  for  X  (  X)  • 


«  S,; 


he 
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The  maximum  energy  output  is  determined  by  the  Planck 


function  (12) 


T  =  the  objects  temperature  in  degrees  Kelvin, 


» 


$ 


h  =  Planck's  constant, 
c  =  the  speed  of  light, 

\=  the  wavelength,  and 
k  =  Boltzmann's  constant  =  1.38E-23 
Joules/molecule-degrees  K. 

Given  a  standard  target,  one  with  a  temperature  of  300 
degrees  Kelvin  and  an  emissivity-area  equal  to  1  square 
meter,  and  a  spectral  irradiance  at  a  distance  r  the  Planck 
function  yields  a  spectral  radiance  L  for  wavelengths 
9  <  X  <  13  microns. 


L,  =  30  Watts-m_2-ster 


■■^■-um  * 


30 


Watts-m“^-um-1 


IT  r2 


For  a  four  micron  bandwidth  from  9  to  13  microns,  the 
irradiance  E  integrated  over  the  spectral  bandpass. 


120 

=  - -  w-m 


-m-2 


T  TT  r^ 

Therefore, 

a.j  -  uo  e  Sn  7aq) 

Kc/i  v  r2 
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Solving  this  equation  for  the  range  r. 


r  -- 


/  u»s  ^a)za)\ 
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The  background  contribution  to  the  signal-to-noise 
equation  must  be  determined.  The  background  current  when 
evaluated. 


e  s«  On _ y(i)TJl) 

hc4 


where  =  the  sensor  wavelength  bandpass. 

This  assumes  7J  (\)  is  constant  across  the  bandwidth  and 

is  the  average  wavelength  (11  microns).  From  calculations 

in  the  previous  section,  the  amount  of  background  radiance 

Lq  at  eleven  microns  varies  from  about  1.0E-11  to  1.0E-13 
“a 

W-cm_^-um-^-ster-^ . 

The  solid  angle  QpoV  is  that  portion  of  the  field 
of  view  observed  by  the  sensor  at  any  instant  in  time. 
Since  the  reticle  slit  only  allows  about  11%  of  the  total 
field  of  view  to  reach  the  detector,  =  0.11  Q  FOy 

For  various  values  of  LB  the  above  equation  gives  the 


values  for  iB  in  table  V. 


TABLE  V 


Sensor  Current 


BACKGROUND  SPECTRAL  RADIANCE  =  La 


CURRENT  ig 


1.0E-11  Watts-cm 5-um  -ster'1  7.37E-11  amps 

1.0E-12  7.37E-12 

1.0E-13  7.37E-13 


The  dark  current  contribution  to  the  signal-to-noise 
equation  must  be  determined.  The  dark  current  is  the 
combination  of  two  factors.  First  there  is  the  detector 
output  current  when  no  photons  are  incident.  Second,  there 
is  the  detector  output  current  from  infrared  radiation 
given  off  by  the  rest  of  the  sensor,  i.e.  optics,  tube, 
baffles,  etc.  By  cooling  both  the  detector  and  the  sensor 
components,  this  noise  can  be  kept  to  a  minimum,  but  not 
entirely  eliminated.  For  the  IRAS  sensor  the  system  in- 
band  noise  equivalent  flux  density  was  1.0E-19  watts-cm-^ 
(20:69).  This  value  included  the  system  noise  from  the 
detectors,  amplifiers,  and  filters.  The  proposed  sensor  is 
about  1/10  a  large  as  IRAS,  therefore,  the  proposed 
sensor's  NEFD  will  be  about  10  times  that  of  IRAS,  1.0E-18 
watts-cm  .  Since  this  value  is  much  less  than  background 
irradiances,  the  dark  current  can  be  neglected  for  the 
accuracy  desired  in  this  investigation. 


m 

$5 


Now  that  the  values  needed  for  the  range  calculation 


are  available,  the  maximum  range  an  ASAT  can  be  from  its 


target  and  still  detect  the  target  can  be  found. 


Assuming  a  signal-to-noise  ratio  of  five,  a  sample 


time  Ts  of  1.39E-3  seconds,  and  minimum  background  current 


of  7.37E-13  amps,  the  equation 


gives  a  value  of  if  =  4.74E-14  amps. 

Solving  for  the  range  equation  below  using  the  values 


from  Table  IV, 


noS,Va)Z(i) 


nV-  rr  iT/ 


Gives  the  maximum  range  of  7271  kilometers.  Table  VI 


lists  ranges  for  other  values  of  background  radiances. 


W" 


TABLE  VI 


Maximum  Ranges  versus  Background  Radiances 
Background  -  Watts-cm^-um'-ster  '  Range-Kilometers 


Spectral  Radiance  r 

IE-13  7271 
2E-13  6143 
3E-13  5563 
4E-13  5184 
5E-13  4907 
6E-13  4691 
7E-13  4516 
8E-13  4369 
9E-13  4244 
IE-12  4135 
2E-12  3482 
3E-12  3149 
4E-12  2931 
5E-12  2773 
6E-12  2650 
7E-12  2250 
8E-12  2467 
9E-12  2396 
IE-11  2333 
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The  values  of  the  maximum  sensor  detection  ranges 
derived  the  the  previous  chapter  depend  upon  the  accuracy 
of  the  Walker  Sky  Model.  While  this  model  is  based  upon 
actual  observational  data,  the  model's  accuracy  may  be 
questioned  from  several  different  angles.  Specifically 
three  questions  are  raised: 

1)  Does  the  model  include  all  the  stars  that  have  a 
significant  contribution  to  the  sky  radiance? 

2)  How  does  the  number  of  infrared  sources  predicted  by 
the  model  compare  with  cataloged  observations? 

3)  What  non-stellar  infrared  sources  must  be  added  to 
obtain  true  sky  radiance  values? 

This  chapter  will  answer  the  first  two  questions.  The  next 
chapter  will  answer  the  last  question. 

Limiting  Ir radiance  Effects 

The  sky  model  data  only  lists  stars  down  to  an 
irradiance  level  of  1.0E-20  Watts-cm~^-um-^.  To  determine 
if  fainter  stars  contribute  to  the  sky  radiance,  the  data 
in  Table  II  can  be  examined.  For  each  galactic  position, 
the  number  of  sources  increases  with  decreasing  spectral 
irradiance.  The  number  of  sources  for  each  entry  is  the 
total  number  of  sources  at  that  particular  irradiance  level 
and  brighter.  It  is  the  number  of  sources  times  the 


irradiance  level  that  determine  the  individual  irradiance 


level  contributions  to  the  total  sky  radiance. 

In  examining  the  columns  of  data,  the  number  of 
sources  increase  by  an  order  of  magnitude  or  more  for  each 
decreasing  irradiance  level  only  at  low  galactic  latitudes. 
At  higher  latitudes,  the  number  of  sources  increases 
rapidly  at  first  and  then  levels  off  as  irradiance  level  of 
1.0E-18  watts-cm“2-um_1  is  reached.  Evaluating  this  trend 
leads  to  the  conclusion  that  the  number  of  sources  appears 
to  be  rather  finite  at  high  galactic  latitudes  and  that 
fainter  irradiance  levels  than  1.0E-20  will  not  contribute 
to  the  model  for  that  sky  region.  Figure  6  shows  this 
trend  by  plotting  number  versus  irradiance.  Figure  7 
shows  this  by  plotting  the  contribution  to  the  total 
radiance  by  each  irradiance  range. 

At  galactic  latitudes  of  zero  and  two  degrees,  the 
number  of  sources  does  not  level  off.  This  indicates  that 
at  fainter  irradiance  levels  than  1.0E-20  watts-cm-2-um~^, 
additional  sources  may  be  of  sufficient  number  to 
contribute  significantly  to  the  model.  For  example,  at  an 
irradiance  level  of  1.0E-21  an  estimate  of  the  number  of 
sources  of  1.0E6  is  possible  if  the  trend  from  irradiance 
levels  1.0E-19  and  1.0E-20  is  followed.  This  is  shown 
graphically  in  figures  8  and  9. 

How  faint  this  analysis  should  be  taken  is  limited  to 
those  irradiance  levels  that  make  significant  contributions 
to  the  model.  The  contribution  for  irradiance  level  1.0E- 


1.0E-20  -19  -18  -17  -16  -15  -14 


Figure  6 
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Irradiance  Watts-cm  -um 


Number  of  Sources  vs  Irradiance  Range 
Galactic  Latitude  20  Longitude  90 


Figure  7.  Change  in  Radiance  vs  Irradiance  Range 
Galactic  Latitude  20  Longitude  90 
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21,  for  latitude  0,  longitude  0,  is  estimated  to  be  as 
follows: 

(1.0E6  -  1.0E5)  x  1.0E-21  =  9.0E-16  W-cm-^-um-^ 

The  radiance  value  for  this  area  was  previous  calculated  to 
be  5.0E-15.  This  means  that  the  estimated  contribution 
from  irradiance  level  1.0E-21  is  about  a  20  percent 
increase  over  the  previously  calculated  value. 

This  process  can  be  extended  to  irradiance  level  1.0E- 

22,  but  this  level  will  require  the  number  of  sources  to  be 
on  the  order  of  1.0E7  if  it  is  to  make  the  same 
contribution  to  the  radiance.  Since  the  Milky  Way  galaxy 
has  an  estimated  hundred  billion  stars  in  it  (1),  expecting 
ten  million  stars  to  be  visible  while  looking  through  the 
galactic  center  is  not  unreasonable.  However,  the  actual 
number  of  stars  detectable  will  be  limited  by  interstellar 
dust  and  gas  clouds.  So  it  is  very  difficult  to  predict 
when  there  are  no  more  fainter  stars  that  will  contribute 
to  the  model. 

In  summary,  the  model  seems  accurate  at  latitudes 
greater  than  two  degrees.  On  the  galactic  plane  the 
radiance  values  previous  calculated  in  Table  IV  could  be  in 
error.  Determination  of  the  magnitude  of  the  error  will 
be  left  to  an  examination  of  actual  observations. 


The  comparison  to  actual  infrared  observations  at 
around  10  microns  is  divided  into  two  parts:  1)  a 


»  •  ,>■  Y 


V” 


comparison  with  data  from  the  AFGL  catalog  data,  and  2)  a 
comparison  with  IRAS  preliminary  results. 

The  AFGL  catalog  of  infrared  stars  has  some  gaps  in 
its  coverage.  These  gaps  are  not  large  at  10  microns,  but 
must  be  considered  when  choosing  areas  of  the  sky  for  the 
comparison.  Four  areas  of  the  sky  were  chosen.  These 
areas  are  15  degrees  on  a  side  and  represent  two  areas 
along  the  galactic  plane,  and  two  areas  at  high  galactic 
latitude.  The  approximate  centers  of  the  areas  are  listed 
in  Table  VII.  The  areas  are  plotted  in  galactic 


coordinates 

in  figure 

10. 

TABLE  VII 

Centers 

of  Selected 

Sky  Areas 

Area 

Galactic 

Latitude 

Galactic 

Longitude 

Right 

Ascension 

Declination 

1 

0 

30 

1900 

0 

2 

5 

210 

0700 

0 

3 

60 

270 

1200 

0 

4 

-60 

90 

0000 

0 

Appendix  B  contains  the  list  of  stars  for  each  area 
with  their  AFGL  catalog  number,  right  ascension, 
declination,  and  magnitude.  Also  in  Appendix  B  is  a 
graphic  plot  of  the  stars  in  each  area.  A  magnitude  of  0.0 
equates  to  8.3E-17  Watts-cm“^-um”^  at  11  microns.  A 
discussion  of  magnitudes  and  irradiance  used  in  this  thesis 


Selected  Sky  Areas  in  Galactic  Coordinates 


is  in  Appendix  C. 

Table  VIII  shows  the  radiance  value  calculated  for 
each  area  based  upon  the  number  of  stars  and  their 
respective  irradiances.  Also  in  this  table  is  the  radiance 
value  predicted  by  the  Walker  Sky  Model.  There  is  general 
agreement  between  prediction  and  AFGL  catalog  data  away 
from  the  galactic  plane,  but  on  the  plane  the  model 
predicts  more  radiance  than  the  catalog  data  indicates. 

This  is  because  the  AFGL  catalog  only  contains  sources  at 
faint  as  1.0E-16  Watts-cm"2-um“1  and  lacks  the  vast  number 
of  faint  sources  in  the  plane  of  the  galaxy. 

It  is  also  possible  to  relate  the  number  of  sources 
predicted  with  AFGL  catalog  data  at  those  irradiance  levels 
for  which  the  catalog  has  data.  The  AFGL  catalog  has  data 
for  irradiance  levels  of  1.0E-14,  1.0E-15,  and  1.0E-15 
Watts-cm-2-um-1.  For  the  four  areas  previously  mentioned. 
Table  IX  shows  the  comparison  between  the  number  of  source 
predicted  from  the  Walker  Sky  Model  versus  the  actual 
number  of  sources  per  square  degree  from  the  AFGL  catalog. 
The  values  of  actual  and  predicted  usually  agree  with  a 
factor  of  ten.  The  best  agreement  occurs  at  an  irradiance 
level  of  1.0E-16  Watts-cm-2-um“*. 

To  summarize,  the  Walker  Sky  Model  does  agree  roughly 
with  11  micron  data  present  in  the  AFGL  catalog.  This  is 
not  surprising  since  the  model  is  based  partly  on  the  data 
in  this  very  catalog.  However,  the  number  of  sources 


predicted  is  an  average  based  upon  data  for  the  entire  sky, 
and  cannot  be  expected  to  agree  precisely  for  any  one 
region. 


TABLE  VIII 

Radiance  Comparison  Measured  in  Watts-cm'"2-um“^-ster“l 
Area  Radiance  Predicted  by  Radiance  from 

Walker  Sky  Model  AFGL  Catalog  Data 


2.1E-11 

1.7E-12 

9.4E-14 

9.4E-14 


2.3E-12 

3.3E-13 

1.2E-13 

7.7E-14 


TABLE  IX 

Sources  per  Square  Degree 
Comparison  of  Predicted  versus  AFGL  Catalog 
Irradiance  Levels  are  Watts-cm-2-um-1 


Irradiance  Level 

Predicted 

AFGL  Cal 

Area  1 

1.0E-14 

5.6E-4 

8.9E-3 

1.0E-15 

1.9E-2 

6.2E-2 

1.0E-16 

7.1E-1 

5.0E-1 

Area  2 

1.0E-14 

5.4E-4 

0.0 

1.0E-15 

1.4E-2 

8.8E-3 

1.0E-16 

1.4E-1 

1.4E-1 

Area  3 

1.0E-14 

1.3E-4 

0.0 

1.0E-15 

2.8E-3 

0.0 

1.0E-1 6 

7.0E-2 

6.2E-2 

Area  4 

1.0E-14 

1.3E-4 

0.0 

1.0E-15 

2.8E-3 

0.0 

1.0E-16 

7 .OE-2 

4.4E-2 

The  comparison  of  the  Sky  Model  results  with  IRAS 
data  is  more  difficult  because  preliminary  data  from  IRAS 
is  all  that  is  currently  available.  This  data  consists  of 
a  small  catalog  of  stars  detected  by  the  IRAS  "minisurvey" 
conducted  before  the  general  observation  program  began 
(22),  and  some  initial  analysis  of  this  data  (29). 

The  IRAS  minisurvey  consisted  of  a  survey  of  300 
square  degrees  between  galactic  latitude  20  to  40  degrees. 
The  data  set  is  considered  95%  complete  down  to  an 
irradiance  level  Of  2.1E-18  Watts-cm“2-um“1  (22:C1).  For 
the  minisurvey  region  at  12  microns,  the  density  of  sources 
brighter  than  2.1E-18  Watts-cm~2-um--*-  is  estimated  to  be 
1.1  per  square  degree  (29:L7).  This  agrees  within  a  factor 
of  two  with  the  value  predicted  by  the  Sky  Model  (see  table 
II)  for  an  irradiance  level  of  1.0E-18  Watts-cm"2-l-um-1 
for  a  galactic  latitude  of  30  degrees  at  11  microns. 
Analysis  of  the  number  of  sources  versus  irradiance  shows 
there  are  more  sources  at  lower  values  of  irradiance  than 
higher  values  of  irradiance,  which  agrees  with  the  sky 
model  (29:L8). 

An  additional  point  of  reference  is  two  scans 
conducted  by  IRAS  that  crossed  the  galactic  plane  (29:L7). 
The  two  scans  were  at  galactic  longitudes  144  and  338 
degrees  respectively.  The  number  of  sources  detected  by 
averaging  these  two  scans  is  shown  in  figure  11.  The 
limiting  irradiance  for  this  figure  is  2.1E-18  Watts-cm“2- 
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Figure  11.  Number  of  Sources  per  Square 

Degree.  IRAS  Data  at  Longitudes 
144  and  338.  Note  general 
agreement  between  this  data  and 
Walker  Sky  Model. 


urn-  .  The  number  of  sources  detected  is  similar  to  those 
predicted  by  the  Sky  Model  for  this  irradiance  level  and  a 
longitude  of  144  degrees,  and  latitudes  greater  than  five 
degrees.  However,  the  number  of  sources  detected  in  the 
center  of  the  galactic  plane  is  about  50%  less  than  that 
predicted  by  the  model. 

The  stars  detected  by  the  IRAS  minisurvey  source 
catalog  are  listed  in  Appendix  B.  The  total  radiance  from 
these  stars  for  the  area  of  the  survey  is  1.1E-13  Watts-cm- 
2-um-1-ster-1.  This  is  a  factor  of  five  less  than  the  sky 
model  prediction  for  galactic  latitudes  between  20  and  40 
degrees.  This  may  be  due  to  the  fact  that  the  data  catalog 
only  contains  stars  which  were  detected  by  IRAS  with  an 
irradiance  of  at  least  2.1E-18  Watts-cm-2-um-1  and  fainter 
stars  can  contribute  to  the  total  radiance  for  that  region. 
The  Walker  Sky  Model  minimum  irradiance  level  is  1.0E-20 

_o  _  ■> 

Watts-cm  4-um  . 

A  reliable  estimate  of  the  radiance  of  the  galactic 
plane  from  the  sky  model  or  from  discrete  source  catalogs 
is  difficult  to  determine.  The  analysis  of  whether  the  sky 
model  includes  stars  faint  enough  shows  that  stars  fainter 
than  1.0E-20  Watts-cm-2-um“^  could  make  significant 
contributions  to  the  total  radiance  along  the  galactic 
plane.  However,  discrete  source  observations  fail  to 
confirm  the  model  prediction  of  radiance  in  the  galactic 
plane.  Price  (25)  has  developed  some  low  resolution  maps 


of  the  galactic  plane  where  total  flux  received  by  the 
sensor  is  plotted  and  not  just  discrete  sources.  The  data 
is  from  the  AFGL  sounding  rocket  flights.  These  maps 
indicate  the  galactic  plane  varies  in  radiance  from  4.7E-11 
Watts-cm“2-um-l-ster~l  at  0  degrees  galactic  latitude  to 
7E-12  Watts-cm”2-um“l-ster~l  at  5  degrees  latitude.  These 
values  are  about  the  same  as  those  predicted  by  the  sky 
model . 


VI. 


to  the  Model 


As  stated  in  the  introduction  to  this  thesis,  sources 
of  celestial  infrared  radiation  include  not  only  stars,  but 
planets,  asteroids,  interplanetary  dust,  and  interstellar 
gas.  Each  contributes  to  the  infrared  background.  It  goes 
without  saying  that  the  Sun  is  the  most  intense  source  of 
infrared  radiation,  but  any  sensor  looking  directly  into 
the  sun  would  be  instantly,  and  probably  permanently, 
blinded. 

All  the  planets  of  our  solar  system  (including  the 
earth)  emit  infrared  radiation  in  the  range  around  ten 
microns.  The  inner  planets  are  the  strongest  sources,  with 
irradiances  falling  as  distance  from  the  sun  increases. 

The  point  source  irradiances  near  the  earth  from  a  few 
planets  at  their  mean  earth-planet  distances  are:  Venus: 
6.22E-12  W-cm-2,  Mars:  3.53E-13  W-cm-2,  Jupiter:  7.40E-14 
W-cm-2  and,  Saturn:  2.24E-17  W-cm-2. 

The  planets  themselves  should  not  cause  many  problems 
for  an  ASAT  sensor.  The  positions  of  the  planets  (the 
earth  and  moon  as  well)  can  be  determined  ahead  of  time 
and  the  ASAT  interception  orbit  adjusted  so  the  ASAT  sensor 
never  sees  a  planet.  If  this  isn't  pocsible,  the  ASAT 
sensor  must  be  "smart"  enough  to  know  when  it  should  see  a 
planet  and  then  ignore  it  as  a  target. 

Asteroids  are  solid  objects  orbiting  the  Sun.  Most  of 


them  lie  in  orbits  between  the  orbits  of  Mars  and  Jupiter. 
Asteroid  surface  temperatures  vary  from  about  180  to  280 
degrees  Kelvin  with  a  mean  of  about  256  degrees  (36:3-31). 
While  this  is  a  little  cooler  than  normal  satellite 
temperatures,  asteroids  still  radiate  a  substantial  amount 
of  energy  at  around  ten  microns. 

Asteroids,  like  planets,  are  confined  to  a  portion  of 
the  sky  along  the  ecliptic.  The  ecliptic  is  a  line  around 
the  sky  inclined  23  and  1/2  degrees  with  respect  to  the 
celestial  equator.  A  graph  of  the  ecliptic  plane  in 
galactic  coordinates  is  shown  in  figure  12.  Ecliptic 
coordinates  are  discussed  in  Appendix  A.  The  asteroids  are 
confined  to  a  band  about  thirty  degrees  wide,  centered  on 
the  ecliptic.  Within  this  band  asteroids  will  raise  the 
value  of  the  infrared  sky  radiance. 

Walker  (36:3-33)  has  developed  a  model,  based  upon  the 
data  from  Murdock  (21),  to  estimate  the  number  of  asteroids 
present  per  square  degree  of  sky.  The  model  was  developed 
from  asteroids  observed  in  the  infrared,  and  orbital  data 
on  1800  known  asteroids.  This  model  for  eleven  micron 
radiation  is: 

N  =  4.6E-19  Ex  exp(-0.14  Be) 

where 

N  =  the  number  of  asteroids  per  square  degree  to  a 
limiting  irradiance  E^, 

E^  =  the  limiting  irradiance  in  Watts-cm“^-um_1,  and 


gurc  12.  Ecliptic  in  Galactic  Coord 


Be  =  the  ecliptic  latitude  in  degrees. 

Evaluating  this  model  for  various  ecliptic  latitudes 
yields  the  data  shown  in  Table  X.  Radiance  values 
calculated  from  this  model  in  Watts-cm-^-um“^-ster-^  is 
shown  in  Table  XI. 

TABLE  X 

Number  of  Asteroids  per  Square  Degree 
Irradiance  Ecliptic  Latitude 


0 

5 

10 

15 

20 

IE-14 

4.6E-5 

2.3E-5 

1.1E-5 

5.6E-6 

2.8E-6 

IE-15 

4.6E-4 

2.3E-4 

1.1E-4 

5.6E-5 

2.8E-5 

IE-16 

4.6E-3 

2.3E-3 

1.1E-3 

5.6E-4 

2.8E-4 

IE-17 

4.6E-2 

2.3E-2 

1.1E-2 

5.6E-3 

2.8E-3 

IE-18 

4.6E-1 

2.3E-1 

1.1E-1 

5.6E-2 

2.8E-2 

IE-19 

4.6E0 

2.3E0 

1.1E0 

5.6E-1 

2.8E-1 

IE-20 

4.6E1 

2.3E1 

1.1E1 

5.6E0 

2.8E0 

IE-21 

4.6E2 

2.3E2 

1.1E2 

5.6E1 

2.8E1 

Irradiance  is  in  Watts-cm“^-um 


TABLE  XI 


Asteroids  in 

Watts-cm“^-um“ 

ic  Latitude 

Radiance 

0 

2.2E-14 

5 

1.1E-14 

10 

5.5E-15 

15 

2.7E-15 

20 

1.3E-15 

& 


The  radiance  for  asteroids  is  brightest  along  the 
ecliptic  plane.  The  contribution  is  at  best  only  about  25% 
of  the  stellar  component  predicted  by  the  Walker  model  at 
galactic  latitude  of  60  degrees  and  higher.  While  the 
asteroid  contribution  to  the  background  may  seem  possibly 
important  for  sky  regions  away  from  galactic  plane,  its 
contribution  is  small  compared  to  the  contribution  from 
zodiacal  light  particles. 

The  infrared  radiation  from  zodiacal  light  particles 
is  thought  to  come  from  particles  that  range  in  size  from 
10  to  100  microns.  These  particles  are  evenly  distributed 
in  the  vicinity  of  the  earth's  orbit  (14:183).  The 
particles  are  at  a  temperature  of  about  250  degrees  Kelvin 
(11:L17).  The  amount  of  infrared  radiation  emitted  from 
the  zodiacal  light  particles  varies  as  a  function  of 
angular  distance  from  the  ecliptic  plane,  and  as  a  function 
of  angular  distance  from  the  Sun  (elongation  angle). 
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The  analysis  of  the  radiance  from  zodiacal  light 
particles  is  not  complete.  The  best  data  available  is  from 
Price,  Murdock,  and  Marcotte  (24:765),  and  from  Hauser,  et 


al  (11:L15).  The  observations  by  Price  are  based  upon 
sounding  rocket  flights  for  the  AFGL  Sky  Survey.  The 
Hauser  data  is  based  upon  scans  made  by  IRAS. 

Both  the  Price  and  Hauser  data  analyses  indicate  that 
the  radiance  in  the  ecliptic  plane  is  much  larger  than  the 
radiance  from  the  stellar  background.  At  an  elongation 
angle  of  75  degrees,  Price  puts  the  infrared  radiance  at  11 
microns  at  about  7.0E-11  Watts-cm-2-um“^-ster~^.  Hauser 
indicates  at  12  microns  the  value  is  about  1.0E-10  Watts- 
cm  •‘-urn  A-ster  A  at  the  same  elongation  angle.  The  main 
difference  between  the  two  data  sets  is  the  value  of  the 
radiance  away  from  the  ecliptic  plane.  The  IRAS  data 
indicates  that  the  radiance  at  the  ecliptic  poles  is  about 
at  least  five  times  as  large  as  the  sounding  rocket  data 
for  the  ecliptic  poles.  Which  data  is  correct  is  still 
open  to  debate.  A  sample  of  data  from  both  sources  is 
listed  in  tables  XII  and  XIII.  The  tables  are  not  at  the 
same  elongation  angle.  A  plot  of  radiances  versus 
elongation  angle  is  shown  in  figure  13.  The  data 
for  figure  13  is  from  Price  (24)  and  Hauser  (11). 
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TABLE  XII 


Zodiacal  Emission 

Infrared  Radiance  Based  Upon  Sounding  Rocket  Data (24) 

Elongation  Angle  of  75  Degrees 

Ecliptic  Intensity  at  }1  Microns 

-ster“* 


Latitude 

Watts-cm“^-um 

-30 

1.1E-11 

-20 

1.7E-11 

-10 

4.3E-11 

-  5 

5.7E-11 

0 

6.8E-11 

5 

5.9E-11 

10 

4.3E-11 

20 

1.8E-11 

30 

1.3E-11 

TABLE  XIII 
Zodiacal  Emission 

Infrared  Radiance  Based  Upon  IRAS  Observations (11) 
Elongation  Angle  of  90  degrees 
Ecliptic  Radiance  at  12  Microns 

Watts-cm“^-um~^-ster-^ 


Latitude 


The  most  important  addition  to  the  sky  model  radiance 
predictions  is  the  radiance  from  zodiacal  light  particles. 
These  particles  are  concentrated  along  the  ecliptic  plane 
but  extend  throughout  the  sky.  The  data  from  Price  (24) 
shows  the  ecliptic  plane  to  be  as  bright  as  the  galactic 
center  (1.0E-11  Watts-cm_2-um“^-ster-^) ,  for  an  elongation 
angle  of  75  degrees.  The  radiance  falls  off  at  higher 
ecliptic  latitude,  but  is  still  seems  brighter  than  the 
stellar  background  at  most  ecliptic  latitudes  (24). 

The  data  from  Hauser  (11)  shows  the  ecliptic  plane  to 
be  as  bright  as  the  galactic  center  (1.6E-10  Watts-cm“2-um~ 
1  -^ter-1).  The  radiance  falls  off  at  higher  ecliptic 
latitudes  as  does  the  Price  data,  but  remains  very  intense 
at  the  ecliptic  pole;  (4.6E-11  Watts-cm~2-um'“1~ster''1) . 

The  Price  does  not  include  values  out  to  the 

ecliptic  pc’es,  1  ut  by  extrapolating  from  the  data 
availab '  ■*>  (table  XII),  the  value  is  at  least  1.0E-12  Watts- 
c”  -a,  -i-ster-1.  This  radiance  is  greater  than  the  sky 
model  prediction  for  galactic  latitudes  above  7  degrees. 

If  the  data  from  Hauser  is  used,  the  infrared  sky  at  12 
microns  is  influenced  mainly  from  zodiacal  emission,  and 
only  the  central  portion  of  the  galactic  plane  adds  any 
contribution  to  background  radiance.  In  fact,  Hauser 
states  in  his  data  analysis  "at  12,  25,  and  60  microns,  the 
zodiacal  emission  is  evidently  the  dominant  source  at  this 
resolution  except  near  the  galactic  plane"  (11:L15).  While 


the  analysis  of  IRAS  data  by  Hauser  is  only  preliminary,  it 
is  suggestive  that  zodiacal  emission  at  wavelengths  around 
11  microns  is  brighter  than  that  estimated  by  Price  (22). 

It  is  difficult  to  plot  radiance  values  for  zodiacal 
light  particles  on  an  all  sky  map  since  the  radiance  varies 
with  elongation  angle  from  the  Sun.  As  the  earth  revolves 
around  the  Sun,  different  stellar  backgrounds  pass  through 
the  entire  range  of  elongation  angles,  0  to  180  degrees. 

It  is  perhaps  best  to  plot  a  median  value  of  zodiacal 
emission  based  upon  a  constant  elongation  angle.  While 
this  will  lead  to  excess  radiance  estimates  at  larger 
elongation  angles,  the  data  is  not  complete  enough  to  be 
very  precise,  however.  Accepting  a  constant  elongation 
angle  will  lead  to  low  estimates  in  the  case  of  smaller 
elongation  angles.  This  can  be  resolved  by  the  fact  that 
the  sensor  must  always  have  an  angular  Sun  avoidance  zone. 
That  is,  a  delicate  sensor  must  not  only  avoid  looking 
directly  at  the  Sun,  but  must  also  avoid  looking  close  to 
the  Sun  because  of  reflections  from  the  optics  and 
structure  of  the  sensor.  An  elongation  angle  of  75  degrees 
is  proposed  since  there  is  data  from  both  Price  and  Hauser 
at  this  angle.  This  simplification  will  mean  that  c  larger 
elongation  angles  the  sensor's  range  will  be  greater  than 
predicted,  but  a  smaller  elongation  angles,  the  sensor's 
range  will  be  smaller. 

The  result  of  the  above  analysis  indicates  that  the 
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infrared  sky  at  wavelengths  around  11  microns  is  divided 
into  two  regions:  the  galactic  plane,  dominated  by  sources 
within  the  central  disk  of  the  galaxy;  and  the  remainder  of 
the  sky,  dominated  by  emission  from  zodiacal  light 
particles  which  are  within  in  the  solar  system. 

From  the  data  from  IRAS,  AFGL,  and  zodiacal  emission 
measurements,  it  appears  that  the  infrared  sky  may  have 
regions  as  intense  at  1.6E-10  Watts-cm~2-um~1-ster-1  (11). 
The  lowest  radiance  could  be  as  low  as  1.0E-12  Watts-cra-2- 
um'-^ster-1  (24),  but  the  minimum  could  be  as  high  as 
4.6E-11  Watts-cm-2-um-l-ster"‘*  (11).  With  this  data  the 
maximum  sensor  range  table  (table  VI)  needs  to  be  adjusted 
to  show  maximum  ranges  for  irradiance  levels  between  1.0E- 
10  and  1.0E-12  Watts-cm“2-um“1-ster“1.  These  ranges  are 
shown  in  table  XIV.  Figure  14  shows  the  possible  ranges 
graphically.  The  procedure  to  calculate  the  ranges  is  the 
same  as  in  chapter  IV. 
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Table  XIV 


Revised  Maximum  Ranges 


Radiance  Level 
Watts-cro*-um  -ster 


1.0E-11 

2.0E-11 

3.0E-11 

4.0E-11 

5.0E-11 

6.0E-11 

7.0E-11 

8.0E-11 

9.0E-11 

1.0E-10 


Range 

km 

2332 

1963 

1774 

1651 

1562 

1492 

1436 

1389 

1349 

1314 


•  Aj 
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A  rough  plot  of  the  minimum  probable  infrared 
background  radiance  around  11  microns  is  shown  in  Figure 
15,  plotted  in  galactic  coordinates.  This  ploc  takes  into 
account  the  estimates  for  both  the  galactic  plane  and 
zodiacal  emission  radiances.  The  maximum  probable  infrared 
background  is  about  an  order  of  magnitude  above  the  values 
in  Figure  15.  There  is  a  considerable  difference  between 
Figure  15  and  the  radiance  contour  chart  in  Figure  2. 


Revised  Radiance  Contour  Map 


Several  conclusions  can  be  drawn  from  the  analysis  in 
this  thesis.  First,  the  data  on  the  infrared  sky  at  a 
wavelength  of  11  microns  is  incomplete.  Existing  sky 
models  are  uncertain  enough  to  dramatically  affect 
satellite  sensor  performance.  Secondly,  more  measurements 
of  the  zodiacal  light  are  needed  to  resolve  the 
uncertainties  of  the  contribution  of  this  source  to  the 
infrared  background.  When  the  entire  IRAS  data  catalog 
becomes  available  and  is  analyzed,  many  of  the  data  gaps 
can  be  filled. 

There  are  several  different  factors  that  could  change 
the  maximum  ranges  that  a  sensor  could  be  from  a  target  and 
still  detect  it  reliably.  These  include  the  size  of  the 
optics,  the  size  of  the  target,  the  detector  material,  the 
number  of  detectors,  and  the  signal  processing  algorithm. 
The  example  sensor  in  this  thesis  used  a  silicon-arsenide 
detector,  much  like  the  ones  used  on  IRAS.  Changing  to  a 
material  like  HgCdTe  could  improve  the  sensor's  ability  to 
detect  fainter  objects  (3).  Increasing  the  number  of 
detectors  to  a  two  dimensional  array  such  as  a  charge- 
coupled  device  and  employing  motion  detection  in  the  target 
signal  processing  could  also  improve  the  sensor's  maximum 
range.  Larger  optics  will  increase  system  sensitivity  as 


well . 


An  area  for  further  study  could  be  an  analysis  of  ASAT 
interception  orbits  to  determine  the  precise  path  across 
the  infrared  background  that  a  sensor  would  cover. 


COORDINATE  SYSTEMS 

This  thesis  used  three  celestial  coordinate  systems, 
equatorial,  ecliptic  and  galactic.  Equatorial  coordinates, 
as  the  name  suggests,  is  referenced  to  the  plane  of  the 
earth's  equator.  Ninety  degrees  from  this  plane  is  another 
plane  in  which  the  north  and  south  celestial  poles  lie. 
Angles  in  the  equatorial  plane  is  measured  in  "right 
ascension,"  using  hours,  minutes  and  seconds,  where  360 
degrees  are  divided  into  24  hours,  60  minutes  per  hour  and 
60  seconds  per  minute.  Angles  in  the  north-south  plane  is 
measured  in  "declination,"  using  degrees,  minutes  and 
seconas  (or  decimal  minutes).  Declination  is  positive 
a: ove  toe  equator,  and  negative  below  the  equator,  just  as 
lat  it  ad*  or  the  earth  is  measured. 

t  *  a .  v  e.n  s  :  ■  >n  is  measured  from  a  zero  point  at  the 
vf  :  •  eg.. x  as  it  is  sometimes  called  "the  first 

..  ’  •  :  A:  e.  ."  This  pc  int  lies  along  a  line  of  the 

**  '  -  ‘  .  r  the  ;  lane  of  the  earth's  equator  with  the 

;  »•  <  :  1  j:u  ore.  t  around  the  sun.  This  point,  the 

>  x,  f  :  era:  ns  fixed  with  respect  to  the 

•  •  ;•  ‘  •  :  *  ras  is  lor  this  coordinate  system.  The 

>  ■  '  •  '  *t  :*  es;  :  n  affects  this  point,  but  only  to 

•  • <  vf  :  ver  y  long  time  spans.  The  angle 

jr,  object  ,  measured  in  the 
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equatorial  plane,  in  a  easterly  direction,  is  the  object's 
right  ascension.  In  figure  A-l,  this  is  the  angle  between 
and  R.  The  angle  from  R  to  the  object,  measured  in  the  R- 
Z  plane,  is  the  objects  declination.  In  figure  A-2,  the 
object  X’s  right  ascension  is  the  angle  OL  and  its 
declination  is  the  angle  5  . 

The  plane  containing  the  solar  system  is  called  the 
ecliptic  plane.  When  describing  the  position  of  objects  in 
the  solar  system  its  is  often  more  convenient  to  define 
their  position  with  respect  to  the  ecliptic.  This 
coordinate  system  uses  ecliptic  longitude,  and  ecliptic 
latitude.  The  prime  meridian  for  this  ecliptic  system  is 
the  vernal  equinox. 

When  the  traces  of  the  planes  of  the  celestial  equator 
and  ecliptic  are  drawn  on  the  sky,  their  point  of 
intersection  is  at  the  vernal  equinox.  These  two  planes 
are  inclined  at  an  angle  of  23  degrees  26  minutes.  This 
angle  is  called  the  obliquity  of  the  ecliptic  6 
Longitude  is  measured  in  degrees  in  an  easterly  direction 
from  the  vernal  equinox.  Latitude  is  measure  in  degrees 
above  or  below  the  ecliptic  plane.  In  figure  A-3,  the 
object  X's  ecliptic  latitude  is  the  angle  Q  ,  its 
ecliptic  longitude  is  the  angle  X 

When  there  is  a  need  to  represent  the  orientation  of 
etars  or  other  objects  in  relation  to  the  galaxy  as  a 
whole,  a  galactic  coordinate  system  is  used.  This  system 


Galactic  Longitude 


tan 


-1 


sir.  (d)  -  sm(b)sin<27.4) 
oos(d)sin<a-192.25)cos(27.4; 


♦  33 


where: 


a  =  the  objects  riqht  ascension 
d  =  the  objects  declination 

The  conversion  of  galactic  to  equatorial  coordinate: 
aocsnpl i shed  by  these  equations: 

Dec!  mat  ion 


-  sin"1 


c  o  s  f  b  )  c  o  s  '  2  7 . 4 )  r.  l  n  (  1  -  3  3 } 
♦  s  l  n  '  b  :  s  i  n  f  2  7 . 4 ; 


Fi  jht  Ascension 


c 

os  '  b )  « : o s  ( 

i-33) 

s  i  n  '■ 

h  )  c  o  s  '  2  7 . 4  , 

♦  19  2.2b 

~  C  O  0  tj 

s  l  n  (  2  7 . 4  / ::  i  n  f  1  -  3  3 

Th<  *•  , river  r;  i  or.  of  equatorial  coordinates  to  ecliptic 
rdmaten  is  accomplished  by  the  following  (6), 


lptic  Lon_.it. do  =  tan 


t  a  n  ( a )  c  o  s  (  £  )  ♦ 


t  a  n  (  d )  s  i  n  f  £  ) 


cos ( a ) 
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Ecliptic  Latitude  =  sin  |sin(d)cos(€  ) 

1  1 

-  cos  (d)  sin  (  £  )  sin(a) 


where 

a  =  the  object's  right  ascension, 
d  =  the  object's  declination,  and 
£  =  23  degrees  26  minutes. 

The  conversion  from  ecliptic  to  equatorial  coordinates 
is  accomplished  by  the  following  (6), 

tan  ( n )  sin  (  £  ) 


Right  Ascension  =  tan 


tan  (  X  )  cos  (  £  )  - 


cos 


(X) 


Decl ination 
where 


X 

p 

€ 


=  sin1  |sin(^J)cos(£  )  +  cos(^])sin(  £)sin  (X)] 

the  object's  ecliptic  longitude, 
the  object's  ecliptic  latitude,  and 
23  degrees  26  minutes. 


Cataloged  Observations 

The  following  data  is  presented  to  two  sections:  the 
cataloged  data  from  the  AFGL  catalog  at  11  microns  (25), 
and  the  discrete  source  data  from  the  IRAS  minisurvey  at  12 
microns  (22).  The  AFGL  data  is  for  four  areas,  each  225 
square  degrees.  The  IRAS  data  covers  300  square  degrees. 
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STAR  DATA  FOR  AREA  1 


Approximate  Center 


Galactic 

Latitude 

0 

Right  Ascension 

1900 

Galact  ic 

Long i t  ude 

30 

Decl inat ion 

0 

AFGL  Number 

Right  Ascension 

Decl ination 

Magnitud 

2210 

183533 

-0650.7 

-6.10 

2271 

185050 

+0110.4 

-5.20 

2390 

192430 

+1115.6 

-4.20 

2349 

191240 

-0708.3 

-4.00 

2381 

192154 

+1424.5 

-3.60 

2223 

183053 

-0221.1 

-3.50 

2206 

183452 

+1024.1 

-3.40 

2368 

191736 

-0806.1 

-3.20 

2402 

192806 

-0254.1 

-3.20 

2461 

194724 

-0743.4 

-3.00 

2177 

182847 

-0207.6 

-2.90 

2251 

184502 

-0203.0 

-2.90 

2334 

190754 

+0900.8 

-2.70 

2749 

195456 

-0200.4 

-2.70 

2139 

182025 

-1342.9 

-2.60 

2290 

185604 

+0638.3 

-2.60 

2350 

191328 

+0934.1 

-2.50 

2376 

192009 

+1358.5 

-2.50 

2227 

183846 

-0424.2 

-2.40 

2241 

184115 

+1353.1 

-2.40 

2324 

190405 

+0807.8 

-2.40 

2341 

191102 

+1047.5 

-2.40 

2147 

182208 

-1316.1 

-2.30 

2169 

182630 

-1055.2 

-2.30 

2178 

182850 

-0838.2 

-2.30 

2213 

183559 

+0845.6 

-2.30 

2310 

190045 

+0724.6 

-2.30 

2067 

180404 

-0941.8 

-2.10 

2113 

181505 

-1146.8 

-2.10 

2165 

182459 

-0351.5 

-2.10 

2190 

183126 

-0720.9 

-2.10 

2245 

184323 

-0242.6 

-2.10 

2291 

185612 

+1256.1 

-2.10 

2379 

192044 

+1410.0 

-2.10 

2119 

181606 

-1357.4 

-2.00 

2132 

181829 

-1304.3 

-2.00 

2304 

185921 

+0107.1 

-2.00 

2345 

191158 

+1104.9 

-2.00 

2088 

180910 

-0435.8 

-1.90 

2117 

181543 

-1346.4 

-1.90 

2154 

182352 

-0655.5 

-1.90 

2256 

184637 

-0658.4 

-1.80 

2272 

185113 

+0036.2 

-1.80 

2343 

191122 

+0003.5 

-1.80 

4250 

193039 

+1337.5 

-1.80 

2440 

194114 

+0338.3 

-1.80 

2136 

181934 

-1331.9 

-1.70 

2222 

183731 

-0236.0 

-1.70 

2223 

183732 

-0545.5 

-1.70 

2284 

185347 

+0751.1 

-1.70 

2378 

192038 

+1423.0 

-1.70 

4237 

182242 

-1318.0 

-1.60 

2316 

190253 

+0809.8 

-1.60 

2333 

190733 

+0920.1 

-1.60 

2374 

191915 

+0923.2 

-1.60 

2398 

192738 

+0249.6 

-1.60 

2152 

182339 

-1151.3 

-1.50 

2157 

182422 

-1242.0 

-1.50 

2205 

183447 

-0527.7 

-1.50 

2314 

190139 

-0546.4 

-1.50 

2108 

181405 

-1211.6 

-1.40 

2150 

182307 

+0543.8 

-1.40 

2203 

183413 

-0738.3 

-1.40 

2380 

192055 

+1447.7 

-1.40 

2423 

193540 

+1138.2 

-1.40 

4252 

194107 

-0004.5 

-1.40 

2057 

180121 

+0826.6 

-1.30 

2192 

183137 

-1133.3 

-1.30 

2199 

183317 

+0532.7 

-1.30 

2202 

183351 

-0723.4 

-1.30 

2243 

184142 

-0423.3 

-1.30 

2400 

192737 

-0056.6 

-1.30 

2065 

180358 

-0814.3 

-1.20 

2083 

180738 

-1033.5 

-1.20 

2127 

181755 

-1348.2 

-1.20 

2186 

183037 

-1410.8 

-1.20 

2193 

183146 

-0845.7 

-1.20 

2207 

183504 

-0622.3 

-1.20 

2266 

184935 

+1207.5 

-1.20 

2287 

185515 

+0322.9 

-1.20 

2288 

185553 

+0435.4 

-1.20 

2371 

191813 

+1349.8 

-1.20 

2161 

182442 

-1200.0 

-1.10 

2174 

182818 

-0945.2 

-1.10 

2185 

183026 

-0730.1 

-1.10 

2211 

183539 

-0532.5 

-1.10 

2259 

184738 

-0758.8 

-1.10 

2276 

185216 

+1035.3 

-1.10 

2392 

192455 

+0656.9 

-1.10 

2453 

194357 

+1030.7 

-1.10 

4253 

194526 

+0921.9 

-1.10 

2118 

181542 

-0655.0 

-1.00 

2168 

182615 

-1134.7 

-1.00 

2188 

183053 

-0910.7 

-1.00 

2194 

183149 

-0759.3 

-1.00 

2229 

183928 

-0505.2 

-1.00 

2070 

180443 

+0633.4 

-.90 

2153 

182352 

-1226.8 

-.90 

2164 

182459 

-0842.6 

-.90 

2166 

182511 

-1304.1 

-.90 

2248 

184426 

-0447.8 

-.90 

2239 

184049 

+1221.7 

-.80 

2326 

190433 

+0705.0 

-.80 

2329 

190540 

+0612.6 

-.80 

2195 

183202 

-0836.1 

-.70 

2142 

182128 

+0335.7 

-.60 

2282 

185341 

-1036.3 

-.60 

2359 

191509 

+1150.9 

-.60 

2123 

181705 

-1220.6 

-.50 

2204 

183444 

-0243.1 

-.50 

2254 

184539 

-0206.6 

-.50 

2179 

182855 

-1000.3 

-.40 

STAR  DATA  FOR  AREA  2 


Approximate  Center 

Galactic  Latitude  5  Right  Ascension  0700 

Galactic  Longitude  210  Declination  0 


AFGL  Number 

Right  Ascension 

Decl ination 

Magnitude 

0877 

060519 

-0623.3 

-2.70 

0915 

061735 

-1036.0 

-2.70 

1028 

065007 

+0827.9 

-2.60 

0865 

060118 
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Approximate  Center 

Galactic  Latitude  60  Right  Ascension  1200 

Galactic  Longitude  270  Declination  0 


AFGL  Number 

Right  Ascension 

Declination 

Magnitude 

1594 

125956 

+0525.9 

-2.50 

1554 

122748 

+0442.8 

-2.20 

1579 

124441 

+0424.8 

-1.70 

1586 

125254 

+0338.6 

-1.50 

1575 

124241 

-0614.9 

-1.40 

1535 

120443 

-0629.0 

-1.30 

1566 

123546 

+0206.2 

-1.20 

1583 

125139 

-0915.8 

-1.10 

1495 

112913 

-1205.3 

-.90 

1549 

122238 

+0101.4 

-.90 

1458 

110453 

-1111.7 

-.80 

1516 

114835 

-1056.2 

-.80 

4156 

123251 

+0618.6 

-.50 

1487 

112320 

+0930.5 

-.30 

rv_ • 

ft 

Iv 


ST  A  5  !  A"  •  3 

A  A  :  A  4 

At  i  :  .  x  ,  -  a  *  c 

*  *  ”  • 

C-i i  act  ic 

La  t  it  .uie  -6  ' 

Gal  act  :c 

Li  r.>j  it  jit  Q  0 

It'.  :  • 

> 

L  Nuru  e : 

F.jht  A.  ae:.;  :  : 

It  .  1*  ■ 

*• 

* 

3  C  9  9 

2  3  2  5  4  ^ 

.  ; ,  3  8.4 

0  C  8  5 

00328" 

-114  6.0 

0066 

002426 

-  C  6  r  4  .  o 

- 

.  4 

3023 

230406 

t  ;  C  ;  t  ' 

.  4 

3  C  3  1 

230658 

+  C  6  2  ?  .  » 

'  • 

30  5  9 

231415 

+1018.3 

- 

3147 

234348 

+0311.3 

- 

' 

3197 

235928 

-0616 . 4 

“  i  . 

3083 

232006 

-1  10"  .  4 

.  "t 

3058 

23141" 

-0  80  1  .  3 

•{$ 

iS 


>>■ 

*6 


86 


Fcl  DATA  FROM  IRAS  MIN  I  SURVEY 
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Appendix  £ 

Infrared  Flux  Units 

There  are  three  flux  units  used  in  this  thesis: 
Watts-cm  ,t-um  ,  magnitudes,  and  Janskys.  Magnitude  is  a 
system  that  originated  with  measuring  stellar  brightness  at 
visible  wavelengths.  The  flux  unit  Jansky  is  used  for 
measurements  of  astronomical  sources  in  the  infrared.  All 
comparisons  in  this  thesis  is  done  in  Watts-cm“2-um“1. 

The  stellar  magnitude  system  was  developed  by  Pogson 
in  1856  (1).  The  basis  of  this  scale  is  the  definition 
that  a  star  of  magnitude  m  is  2.512  (or  the  5th  root  of 
100)  times  as  bright  as  a  star  of  magnitude  m+1.  That  is, 
the  higher  the  magnitude  of  the  fainter  the  star.  For 
negative  magnitudes,  a  more  negative  magnitude  is  brighter 
than  a  less  negative  magnitude.  If  the  light  flux  from  two 
stars  is  Ll  and  L2,  the  magnitude  (abbreviated  m) 
difference  between  them  is: 

/  L2 

ml  -  m2  =  2.5  log  [ - 

VL1 

Table  C-I  shows  magnitude  differences  and  approximate  light 
ratios.  For  magnitudes  listed  in  Appendix  B,  the  reference 
magnitude  of  0.0  is  set  at  8.3E-17  Watts-cm""2-um”*  (27). 

The  Jansky  is  a  flux  unit  defined  as  1.0E-26  Watts- 
meter-2-Hertz-*  (16).  This  unit  is  wavelength  dependent. 
The  Jansky  is  used  as  a  flux  unit  for  the  IRAS  data  in 
Appendix  B.  The  wavelength  for  the  IRAS  data  presented  is 


Appendix  Q. 

Walker  Sky  Model  Derivation 


This  appendix  shows  the  derivation  of  the  formulae  for 
the  Walker  Sky  Model.  The  model  predicts  how  many  infrared 
sources  can  be  seen  in  a  1  by  1  degree  area  of  the  sky  in 
the  direction  defined  by  the  galactic  latitude  b  and 
longitude  1. 

The  number  of  sources  seen  along  any  line-of-sight 
will  be  a  function  of  the  density  of  the  sources.  The 
Walker  Model  assumes  there  are  three  classes  of  objects, 
dwarf  stars,  giant  stars,  and  supergiant  stars.  The  model 
assumes  that  density  of  stars  of  type  i  at  any  particular 
point  (r^,  b)  is  given  by 


where 


n^  =  the  mean  space  density  along  the  galactic 
plane  in  numbers/cm^, 
r^  *  the  distance  from  the  sun, 

Rc  =  the  radius  of  the  galactic  nucleus, 

=  the  distance  in  the  galactic  plane  from  the 
galactic  center  to  some  arbitrary  point 
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defined  by  (see  figure  1) , 


b  =  the  galactic  latitude 
Oi  =  the  standard  deviation  of  the  galactic 
distribution  of  a  particular  class  of 
objects  in  terms  of  perpendicular  distance 
from  the  galactic  plane. 

The  number  of  objects  seen  in  a  1  by  1  degree  area 
will  be  determined  by  the  number  of  objects  in  the  volume 
formed  by  projecting  the  1  by  1  degree  area  to  maximum 
distance  Rmax.  This  is  shown  in  figure  D-l. 

The  differential  volume  element  is 

A V  s  Ar  (Abr)Z 

where 


r  =  the  distance  from  the  sun  to  the  volume 
element, 

A r  =  the  length  of  the  volume  element, 

Ab  -  the  angular  width  of  the  volume  element. 
The  total  number  of  objects  in  the  volume  is  the 
density  of  the  objects  times  the  total  volume, 


f/  71.  f  (Ab)  (L C- 


where 


Rmax  =  the  maximum  distance  that  objects  can 

still  be  detected  by  a  sensor  before  they 
become  too  faint, 


The  distance  Rmax  is  dependent  upon  the  irradiance  of 
the  sources  and  cannot  be  measured  directly.  However,  a 
minimum  detectable  radiance  can  be  defined,  and  a  change  of 
integration  variable  is  in  order.  Since  the  radiance 


then  z 


where  1^  L  =  the  mean  spectral  intensity  of  the  ith 
class  of  objects  in  watts/ster. 


Then  the  limits  of  integration  become 


1 

o 

If 

—  ^  -  - 

r~  fW  _ 

— »  *  K) 
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2 

Substituting  these  terms  for  r^  and  dr^  leads  to 


Since  I.  ■  and  nj  are  constants, 

N--  I  J";/'  FE;%^ 

M.lJ 


The  total  number  of  sources  is  the  sum  of  all  classes, 


Substituting  for  F  gives 


Since  *  L  ■  /E.  and  allowing  k  *  1.52E-4  to  convert 

*■  A/a  a 

from  number/ster  to  nuraber/deg^,  Walker's  formula  is 


arrived  at 
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